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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 

1)  As  Integrated  circuits  are  miniaturized,  the  surface  becomes  a  large  fraction 
of  each  device,  and  Its  Influence  on  device  performance  becomes  greater.  Moreover,  a 
detailed  understanding  of  semiconductor  surfaces  Is  needed  for  the  fabrication  of 
heterostructures  and  superlattice -based  devices.  Because  very  little  Is  known  about 
the  surfaces  of  III-V  compound  semiconductors,  we  have  launched  a  major  program  to 
understand  their  surfaces  and  native  defects  on  those  surfaces.  This  report  focusses 
on  InAs,  InP,  and  InSb,  and  the  defects  we  have  observed  on  their  (110)  surfaces.  To 
our  knowledge,  these  are  the  first  observations  of  such  step- structure  defects  on 
these  surfaces . 

We  have  been  among  the  first  to  obtain  real -space  atomic -resolution  scanning 
tunneling  microscopy  images  of  the  (110)  surfaces  cf  InAs,  InP,  and  InSb  cleaved  in 
ultra-high  vacuum  (UHV) .  We  observed  steps  of  several-  and  one-atom  height  on  these 
surfaces,  as  well  as  mages  of  single  atoms.  Fig.  1  shows  Individual  In  atoms  on  a 
flat  InAs  (110)  surface,  imaged  under  positive  bias  (which  reveals  only  the  In 
atoms).  Each  bump  is  an  In  atom  spaced  from  Its  neighbors  by  several  Angstroms.  In 
addition  to  flat  terraces,  we  also  observed'  steps  In  three  principal  directions; 
along  (110),  (111),  and  (112)  directions.  An  example  on  the  UHV-clcaved  In?(ll0) 
surface,  showing  all  three  types  of  steps  Is  given  In  Hg.  2,  where  the  three  types 
of  steps  meet  to  form  a  raised  triangular  terrace. 

It  Is  widely  believed  (but  not  directly  demonstrated)  that  atomic- layer-by-layer 
growth  by  molecular  beam  epitaxy  (and  other  methods)  is  Initiated  at  such  steps  on 
(001)  surfaces.  Therefore  we  have  tried  to  understand  the  structure  and  dynamics  of 
these  steps.  We  observed  periodic  variations  in  geometric  structure  along  some  of 
these  steps,  suggesting  that  the  step-edges  reconstruct.  Therefore  we  have  proposed 
models  of  the  nature  of  the  steps  and  their  reconstructions  In  these  cases.  Fig.  3 
shows  an  example  of  the  observed  periodicity  along  a  (112) -type  step  on  InAs(llO). 
Our  atomic  model  of  this  step  Is  presented  In  Fig.  4  and  Fig.  5.  The  experimental 
data  Indicate  a  periodicity  of  6.7±0.6A.  In  the  model,  which  proposes  a  onc-atom-hlgh 
step,  dimers  form  at  the  step  edge,  between  the  top  and  the  bottom  terrace.  The 
theoretical  periodicity  Is  7.4  A,  In  agreement  with  the  data.  We  have  also  proposed 
a  similar  model  of  the  (111) -type  step  edge  and  Its  reconstruction.  Thus  we  have 
laid  the  foundation  for  studying  steps  and  their  roles  In  growth  and  surface 
reactions . 

2)  We  are  continuing  our  studies  of  STM  Imaging  with  a  semlcondcutor  tip, 
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rtportad  pr«vio\isly. 

3)  Ue  are  studying  the  effects  of  electrical  current  (and  its  direction)  on  step 
edges  and  the  domain  populations  of  externally  strained  Si (100)  2x1.  These  are 
atomic -level  electromigration  studies.  These  experiments  are  being  done  in 
collaboration  with  Dr.  Ig  Tsong.  He  is  providing  the  graduate  student  and  we  are 
providing  the  equipment. 

4)  In  October  of  1992  Dr.  Ruth  Nicolaides  of  Picatinny  Arsenal,  in  Dover,  New 
Jersey  visited  as  part  of  an  on-gong  collaboration. 

5)  Our  ARO  research  has  benefitted  from  interaction  with  Pachyderm  Scientific 
Industries,  a  spun-off  company  that  deals  in  STM  hardware  and  software. 


Copies  of  abstracts  of  each  paper  are  attached. 

There  has  been  continuing  communication,  including  visits  by  Uilliam  Packard  to 
Dr.  Ruth  Nicolaides  of  Picatinny  Arsenal,  and  by  Dr.  Ruth  Nicolaides  from  Picatinny 
to  our  program. 


Our  File:  report. arol93.rno 


Fig.  1.  Scanning  tunneling  microscopy  image  of  the  fiat  InAs(llO)  surface  cleaved  in 
ultra-high  vacuum.  The  sample  bias  is  -)-0.5  V,  and  so  only  the  In  atoms  are  visible.  The 
tunneling  current  is  100  pA,  and  the  image  size  is  40  A  x  18  A. 
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Fig.  2.  Triangular  terrace  observed  by  scanning  tunneling  microscopy  of  the  InP(llO) 
surface.  This  terrace  is  formed  by  intersection  of  (110),  (111),  and  (112)  steps.  The  image 
is  450  A  X  450  A,  taken  at  a  bias  of  +2.0  V  and  a  current  of  100  pA. 


Fig.  3.  Scanning  tunneling  microscope  image  of  a  one-atom- high  (llO)-type  step  on 
the  InAs(llO)  surface.  The  observed  periodicity  along  the  st^  edge  is  6.7±0.8  A,  and  is 
attributed  to  dimerization. 
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Fig.  4.  Dimerization  model  for  a  (112)-type  step  edge  on  InAs(llO),  a  [112]  step, 
before  any  reconstruction  or  relaxation,  (a)  Top  view  of  the  step,  (b)  Side  view,  showing 
dangling  bonds  as  “rabbit  ears”  prior  to  dimerization. 
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Electronic  contribiitioas  to  tcBBnipg^tnnneUag-microtcopy  images 
of  an  annealed  ^iC(lll)  sniface 


M.>H.  Tsai,  C.  S.  Cban|,  Jobs  D.  Dow,  and  L  S.  T.  Tsong 
Dtpermmt  of  tkytia,  Aritona  Stmt  Umataity,  Ttmpt,  Araoiu  iS3S7‘1304 
(Rccdvtd  13  July  1991) 

Tbc  raoonstruetion  of  ibt  cubic  B-SiOl  1 1)  luifscc  was  studied  after  aancalini  at  — 1200 *C  The  tur* 
iaoc  coosUtently  showed  a  teomeuy  when  measured  by  lew<acrgy  clectroa  difiaetion 

(USED)  aud  a  6X6  leometry  wbeo  itnaged  by  scanoing  tunaelisi  microscopy  (STM).  To  rmolve  the 
diserapaacy,  we  carried  out  fim>priaciples  cakulatioas  of  the  clectrouie  structure  of  a  surface  model  ia 
which  a  graphite  inooolayer  was  iDcommeasuiatcly  grown  on  top  of  the  Si-ienainated  ^SiOl  1 1)  sur* 
face.  The  energy  band  and  dcasity-of-states  diagruas  provide  an  explanation  to  the  observed 

differences  in  LEED  and  STM  results.  They  also  ezplaio  the  voltage-depeadent  contrast^revena) 
pbenotnenon  observed  in  the  STM  images  of  the  annealed  surface. 


L  INTRODUCTION 

Scanaing>tttnneling>microscopy  (STM)  images  are 
known  TO  sample  the  local  denuty  of  elecuooie  states 
near  the  Fermi  level  at  a  surface.  In  eenain  cases,  these 
images  can  be  associated  with  the  atomic  geometry  of  the 
surface,  the  most  well  known  of  which  is  the  Sid  11)* 
7X7,'  provided  that  the  surface  electronie  structure  is 
closely  tied  to  individual  surface  atoms.  In  other  cases, 
image  interpretations  are  not  so  straightforward,  e.g.,  in 
Si(lll)-(vixv^)Al,*  O  on  GaAs  (110),^  and  Si(lll)- 
(V^XV^)Ag.*'*  Nevertheless,  in  the  majority  of  cases, 
STM  images  do  reflect  the  surface  geometry  and  peiiodi* 
city  shown  in  LEED  patterns.  In  the  course  of  our  $tud> 
ies  on  the  reconstnietions  of  the  euhie  ^StCdll)  sur* 
face,**^  however,  we  found  a  disagreement  between  the 
surface  geometry  determined  by  STM  and  that  deter* 
mined  hy  LEED.  STM  images  show  a  6X6  geometrv 
while  the  LEED  patten  indicates  a  6'OX6v^ 
geometry.  To  explain  this  discrepancy,  we  have  invoked 
a  model  originally  proposed  by  van  Bomme],  Crombeen, 
and  Van  Tooren.*  Three  successive  carbon  layers  of  ^ 
Sic,  after  evaporation  of  silicon  hy  annealing  to  1200*C, 
collapse  into  one  single  laver  of  carbon  atoms  with  a  sur* 
face  density  of  3.67  X 10'^  atoms  em''^  very  elose  to  the 
density  of  a  graphite  monolayer  of  3.80X10'^ 
atoms  cm‘~^.  This  layer  of  carbon  atoms  then  have  only 
to  be  displaced  over  sauE  distances  to  form  a  graphite 
layer.  The  6v^  X  6v^  LEED  patten  is  due  to  multiple 
seatterihg  from  the  graphite  monolayer  on  top  of  the  Si* 
terminated  ^-SiCd  1 1 )  surfaee.^'^ 

In  this  paper,  we  present  first-principles  caleulations  of 
the  electronic  structure  eff  the  annealed  ^iC(lll)  sur* 
face  twsed  on  the  above  mode]  in  order  to  explain  the 
6X6  geometry  observed  in  STM  images.  Because  of  the 
large  number  of  atoms  involved  in  a  unit  ceB,  calcula¬ 
tions  were  carried  out  for  the  graphite  monolayer  and  the 
Si-terminated  ^SiQUl)  surface  separately.  A  compar¬ 
ison  of  the  density  of  sutes  of  these  two  surfaces  is  shown 
to  produce  the  6  X  6  geometry. 
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n.  EXPERIMENT 

A  description  of  our  STM  with  iu  associated  scanning 
tunneling  spectroscopy  (STS)  system  has  been  given  pre* 
viously.***  The  (111)  surfaces  of  euhie  B-SiC  monocrystal¬ 
line  films  were  grown  on  the  Si-terminated  (0001)  siirfaee 
of  hexagonal  a*SiC  suhstrates  hy  the  technique  of  chemi¬ 
cal  vapor  deposition  (CVD).’  Prior  to  tunneling,  the 
sample  surface  was  annealed  to  —llOO'C  for  1  min  to 
eliminate  the  surface  oxide  layer.  After  this  treatment, 
the  surface  consistently  showed  a  6v^  X6v^  LEED  pat¬ 
tern,  in  agreement  with  the  results  of  van  Bommel, 
Crombeen.  and  Van  Tooren.*  All  STM  images  of  the  an¬ 
nealed  surface  were  acquired  at  a  pressure  of  5  1 X 10''** 
Torr  at  room  temperature. 

m.  RESULTS 

Figure  1  shows  voluge-dependent  scans  of  the  same  re- 
gioo  on  the  annealed  B-SiCdll)  surface.  The  surface 
shows  a  6X6  geometry  with  a  unit-cell  edge  of  19  A, 
which  is  approximately  six  times  the  unit<ell  edge  of  3. 1 
A  of  the  ^SiCd  1  IHl  X 1 )  hulk-terminated  surface.  The 
images  (a)  and  (b)  of  Fig.  1  were  taken  simultaneously  on 
the  same  area  with  the  tip  biased  at  3.7  V  when  scan¬ 
ning  from  left  to  right  and  biased  at  —3.7  V  when  scan¬ 
ning  from  right  to  left.  The  hysteresis  in  the  scans  was 
small  and  had  been  correaed  in  the  images  shown  here. 
The  most  remarkable  feature  in  these  two  scans  is  the  oc* 
eurrenee  of  a  contrast  reversal.  The  honeycomb  panem 
in  the  filled'SUie  image  (a)  turu  bto  a  centered-hexagon 
array  in  the  empty-state  image  (bk  The  miniTna  of  the 
honeycombs  in  (a)  comcide  exactly,  as  far  as  the  experi¬ 
ment  can  tell,  whh  the  maxima  of  the  centered  hexagons 
in  (b).  Each  honeycomh  minimuin  and  eaeh  hexagon 
maximum  has  a  diameter  of  — 12  A,  clearly  mueh  larger 
than  the  sixe  of  a  single  atom. 

Figures  2(a)  and  2(b)  show  filJed-state  and  empty*$ute 
images  taken  with  tip  biases  of  *r3  and  —3  V,  respective¬ 
ly,  on  the  same  area.  While  the  eorrugations  of  filleo- 
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FIG.  r  STM  iffiafcs  of  ihc  ^-SiCOll.*  surface  afxcr  anneaici  at  icmperiiurc  ^120C*C*  Scar*  area  is  fOCXSOC  A".  fVcnically*  1- 
di'  =  a:  honzonully*  \  div*  ^00  Aj  TuTii^chnf  curfeu;  ^  <  nA.  Uf  FiViea^suie  image  a**  ;;p  bias  F.  V,  anc  ib*  c?rip*A*suit 
imagr  a;  up  bias  r.  =  —  j.T  V. 


s;a;t  image  la'  appear  unchanged  compared  to  iht  terre* 
sponcinc  image  in  Fie.  lu.-.  the  corrugations  in  the 
empty-state  image  (b;  have  diminished  considerabiy.  In 
Fig.  I*,  we  plot  the  corrugation  as  a  function  of  tip  bias 
for  mages  uker.  ai  different  biases.  1:  is  obvious  that  the 
corrugations  at  pcsitjvc  tip  bias,  i.e.*  filled-sutc  images, 
remain  rciaiivcly  consuni,  white  those  at  negative  bias, 
i.e..  tmpty-sraic  images,  decrease  as  the  tip  bias  de* 
creases.  In  lac:,  a;  tip  bias  less  than  —2.5  V.  the  eor.^u^ 
gations  aimost  compieiely  disappear. 

rv.  Disci-ssios 

We  have  previously  reponee*  the  use  of  STM  and  STS 
tc  show  ihai.  upon  annealing  lo  1200 ’C,  a  graphite  layer 
IS  incommensuraiely  grown  on  top  of  the  Si- terminated 
^-SiC.Tli*  surface.  Our  resui:  agrees  with  previous  re¬ 
sults  obtained  by  other  workers  using  LEED.  .Auge.'  cicc- 
tion  spectroscopy.  ^-lay  phoioemission  spectroscopy. 


and  cifc:ron-enefgy-loss  spectrcsccpy**^^'”*"  on  he.xcgc- 
nal  c*SiC(000n  surfaces  which  are.  of  ccu.^se.  equiv^ien; 
to  the  cubic  iS-SiC(l  3 1}  suriace.  In  Fig.  4,  we  have  drawn 
a  g*^a5hitc  mesh  wmh  a  penocic:r>*  tor  unit-cell  edge  I  o: 
1.46  A  over  a  Si  bulk-tenmnatec  11X1)  surface  of  iS- 
SiGlI  1)  with  a  periodicity  of  i.!  A.  Toe  positions  where 
the  top  C  atoms  coincide  exactly  with  the  Si  atoms  un¬ 
derneath  are  marked  with  circles  and  labeled  A,  £,  C 
£*  and  F.  These  positions  give  nsc  to  the 
geometry*  observed  by  LEED  since  the  unit-cell  edges. 
AE  and  C£*  each  have  a  length  of  13  graphite  hexagons, 
i.e..^  13x2.46  a*  which  is  annually  identical  to 
6^3  X  3.1  A.  Around  each  of  the  positions  A-f.  the.-c 
arc  three  pairs  of  C  and  S:  atoms  at  apexes  of  the  dotied 
tnangic  lying  very  cbsc  to  each  othc.',  i.e.  wiihir  C.3  .A 
laicrally.  The  position  G,  giver,  by  the  donee  circie  in 
Fsg.  4,  ooes  not  have  any  C  atoms  directly  abovt  a  S: 
aioir*.  Six  C  aioms  on  the  doitcc  circie.  howeve.n  arc  n- 
ing  very  close  to  six  Si  atoms,  speculated'  that 


FIC-  1.  STM  images  of  the  S-SiC.:!!:  sunacc  after  annealed  at  icmorniure  --1200‘C.  Scan  area  is  SOOXSOC  A'.  Tcr.nrlmg 
currr.v.  :s  1  .-A.  la  Friicd-staic  imacf  a;  up  bias  1*  ■*>  V*ind  ib‘  enpi) -stair  image  at  tip  bias  F.  »  —  V. 
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ABCF  forms  the  6x6  gcomcin*  ohscr\*cd  by  STM.  To 
expUio  the  contrast  reversal  phenomenon  as  the  up  htas 
is  reversed,  we  need  detailed  knowledge  of  the  electronic 
siniciure  of  the  surface. 

To  calculate  the  6X6  surface  directly  is  beyond  the  ca¬ 
pability  of  a  firsi-princtplcs  method  because  each  uni;  cell 
contains  hundreds  of  atoms.  Thus,  our  approach  to  solve 
the  problem  ts  to  model  the  system  separately  hy  a  graph¬ 
ite  monolayer  and  a  Si-termtnated  ^-SiC(lll)  surface. 


FIO.  4.  A  model  of  a  graphite  monolayer  ^represented  by 
honeycombs)  incommcnsurately  grown  on  the  (1X1)  Si. 
terminated  /^-SiOlll)  surface  (represented  by  dou).  The 
dasbed  mangles,  denoted  A,  B,  C.  D,  £.  and  F.  indicate  posi¬ 
tions  where  the  graphite  atoms  at  the  centers  of  mangles  coin¬ 
cide  exactly  with  the  Si  atoms  in  the  second  layer,  and  those 
graphite  atoms  at  the  apexes  almost  coincide  with  the  Si  atoms 
underneath.  Tne  dashed  circle.  C.  outlines  positions  adhere  the 
graphite  and  Si  atoms  nearly  coincide.  The  unit-cell  spacing  is 
2.46  A  for  graphite  and  j.l  A  for  the  (!  X  1 1  Si-tcmunaicc  sur¬ 
face. 


Wc  use  the  firsi-principlcs  pseud  ofuncuor.  (PSF) 
method in  a  local-denstty  approximation  (LOA)  (Ref. 
14)  to  calculate  the  electronic  structures  of  the  graphite 
monolayer  and  the  ^StCdl  1)  surface. 

For  the  ^SiCClll]  surface,  an  isolated  eight  atomic* 
layer  (or  four  douhle*layer)  slab  model  is  used.  Because 
of  the  double-layer  nature  of  ^-SiCdll),  one  surface  of 
the  crystal  ts  terminated  in  Si  and  the  opposite  surface  is 
terminated  in  C  Each  C  atom  on  the  opposite  surface  of 
the  slab  is  attached  with  a  hydrogen  atom  to  saturate  the 
dangling  bond  in  order  to  prevent  charges  from  transfer¬ 
ring  hack  and  forth  between  the  two  surfaces  of  the  slab, 
thus  making  convergence  possible.  ^Wc  use  the  experi¬ 
mental  lattice  parameter  of  4.3596  and  the  hulklike 
crystal  structure  for  the  5-SiC.  For  the  graphite  mono¬ 
layer,  we  also  use  an  isolated  slab  model.  Tne  lattice  pa¬ 
rameter  is  2.436  A  as  given  in  the  Iheraturc.***^’  The  spe¬ 
cial  k-potnt  scheme  of  Cunningham^*  for  a  two- 
dimensional  hexagonal  lattice  is  used  to  approximate  the 
integration  over  the  hrst  Brillouin  zone  in  calculating  the 
self-consistent  charge  densities  and  potential  for  both 
cases,  ^  t  use  one  spectal  k  pom;  for  ^SiC  anc  six  spe¬ 
cial  k  points  for  graphite. 

Calculated  energy  hands  for  the  f-SiC(lll)  and  the 
graphite  monolayer  are  shown  in  Figs.  3  and  6,  respec¬ 
tively.  In  both  figures,  energies  are  given  relative  to  the 
vacuum  level.  In  Fig.  3,  there  is  a  very  flat  half-fiUcc  Si 
dangling-bond  energy  hand  around  -^3  eV  inside  the  fun- 
damenud  energ>*  gap,  giving  rise  to  a  very*  strong  peak  in 
the  panial  density  of  .>utes  of  the  Si  surface  atoms  as 
shown  in  Fig.  7.  Tne  calculated  energy  gap  is  !.43  eV, 
which  is  about  half  of  the  experimental  value  2.3  eV.^* 
This  discrepancy  is  weU  known  to  he  the  result  of  a 
deficiency  of  the  LDa.*’  The  Fc  level  lies  within  the 
dangling-bond  sute  hand  and  is  4.S  eV  below  the  vacuum 
level,  i.e.,  the  work  function  is  4.S  cV. 

For  the  graphite  monolayer,  the  energy  hands  shown  in 
Fig.  6  are  similar  to  those  of  Tatar  and  Rabii^^  using  a 
modified  Korringa-Kohn-Rostokcr  approach  for  three- 
dimensional  graphite  and  also  to  those  of  Posiemak 


SiC{l11)  energy  Danes 
Si  terntmaiae 


FIG.  3.  Calculated  enerfy.band  diapam  of  the  Si-termL-iaiec 
surface  of  /J-SiCd  1 11.  Dashed  line  indicates  the  Fc.^i 


Monoliyif  ytp***^ 


FIG.  6.  Caictflited  cDer^-band  diafnm  of  the  graphiit 
mo&olayer. 

er  using  the  full-potential  linearized  augmented- 
plane-ti*ave  fFLAPW)  method  for  a  graphite  monoityer* 
We  obtain  a  work  function  of  4»2  eV  for  the  graphite 
monolayer  compared  to  5.D2  eV  obtained  by  Postemak 
er  o/.^  Previous  calculated  work  functions  scatter  over  a 
range  between  5.0  and  8.0  eV,^^'  The  eapcrimental 
work  function  for  bulk  graphite  is  S.O  eV*^*  The  work 
function  for  the  hulk  material  is  expected  to  be  higher 
than  that  for  the  monolayer  because  the  escaping  electron 
experiences  an  extra  attractive  potential  from  the  inner 
layers  despite  the  fact  that  the  couplings  between  graph¬ 
ite  layers  are  weak.  Thus  our  calculated  work  function 
for  xbt  graphite  monolayer  appears  reasonable. 

Tnt  contributions  to  the  STM  images  can  be  approxi- 
maiely  divided  into  two  parts>  one  from  the  Si  dangling- 
bond  states  and  tbe  other  from  the  graphite  electronic 
sutes.  The  Si  dangling-hond  band  is  half  filled.  Howev¬ 
er,  tbe  band  and  the  corresponding  density  of  sutes  are 
not  s>'mmetrical  with  respect  to  tbe  Fermi  level  as  indi¬ 
cated-in  Fig.  5.  Referring  to  Fig.  5,  the  fiUed-sute  pan, 
i.e.,  the  pan  below  the  Fermi  level  (dashed  line),  Is  ex- 


VaojumiMi 


FIG.  7.  Calculated  panial  density  of  suies  for  the  Si¬ 
lt  munated  surface  of  ^-SICU 1 1 }  and  for  the  graphite  inonola>er 
froc:  Ref.  16.  The  Fermi  levels  with  respect  tp  each  Ptbc:  arc 
shomu  below  tbe  vaeuuzn  level. 


trtZDdy  flat  and  is  in  the  large  wave*vector  legioiL  In 
cootrast,  the  empty-ftate  pan  surrounding  tbe  F  point 
has  a  more  pronounced  curvature,  Lt^  larger  energy 
spread,  and  has  smaller  wave  vectors*  Thus  the  empty 
sutes  arc  more  delocalized  and  the  corresponding  density 
of  sutes  is  much  smaller  than  that  for  ^e  filled  sutes. 
Because  of  this,  one  expects  the  contribution  of  the  Si 
dangling-bond  sutes  to  be  larger  for  fiUed-sute  images 
(i.e*,  positive  tip  bias)  than  for  empty-sute  images  (nega¬ 
tive  tip  biash 

Since  the  work  function  for  the  graphite  monolayer  is 
smaller  than  that  for  the  Si-terminated  B-SiC(Hl),  there 
is  a  net  charge  transfer  from  the  graphite  to  the  B^StC. 
Tbe  graphite  Fermi  level  is  located  at  the  r  band,  which 
is  composed  solely  of  p.  orbitals;  here,  r  is  tbe  coordinate 
normal  to  the  surface.  When  there  is  a  Si  atom  directly 
beneath  or  close  to  a  C  atom  in  the  graphite  layer,  the 
charge  transfer  from  the  C  atom  to  underneath  tbe  Si 
atom  will  be  maximum*  Otherwise,  the  charge  transfer 
will  be  small  due  to  the  highly  directional  p.  orbiul.  W‘e 
write  p\iE)  for  tbe  density  of  sutes  of  those  C  atoms 
which  have  a  St  atom  directly  underneath  or  dose  to 
them,  i.e..  C  atoms  in  regions  A-G  in  Fig.  4.  and  P:(r) 
for  the  density  of  all  other  C  atoms.  Tbe  density  of  sutes 
for  hulk  graphite  has  been  given  by  Tatar  and  Raht.^^ 
Since  the  energy  hands  of  tbe  graphiie  monolayer  are 
similar  to  those  of  bulk  graphite,  the  density*  of  sutes  for 
them  will  be  similar*  W'itbin  about  2.5  eV  of  tbe  Fermi 
level,  the  density  of  sutes  is  approximately  parabolic, 
which  is  shown  in  Fig.  8.  Thus 

.  (1) 

where  X  is  a  proportionality  constant. 

For  those  C  atoms  with  a  Si  atom  underneath  or  in 
close  proximity,  charge  transfer  from  C  to  Si  atoms  will 
shift  the  density-of-sutes  curve  dowTi wards  by  an  amount 
6  as  shown  in  Fig.  8,  resulting  in  an  attractive  elecirosui- 


FIG.  8.  Shift  in  tbe  densiiy  of  states  of  the  graph] it  mono¬ 
layer  (Ref.  1^}  as  a  result  of  charge  transfer.  Tne  number  of 
Stales  samples  by  different  tip  biases  is  also  shown. 
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ic  poiendal.  5  has  a  positive  value  in  cncrf  y  units.  Thto 

p,(£)*^/:(£-£,r  +  6)^  .  (2) 

For  lip  bias  V,,  let  A^j(£/-,  and  ^ 

numbers  of  lUtes  between  £^  and  for  the  kinds 

of  C  atoms  mentioned  earlier.  Let  he  the  difference 
in  the  number  of  states  sampled  by  STM  from  the  two 
kinds  of  C  atoms  on  the  surface,  which  in  turn  fives  rise 
to  the  contrast  observed  in  Fifs.  1  and  2.  Hence 

.  (3) 

We  write  for  F,  >0  or  positive  lip  bias  (sampling 
filled  suics).  and  niv*“  for  F,  <0  or  negative  tip  bias 
(sampling  empty  sutes).  From  Eqs.  (1)-(3K  we  now 
have 

^l'_^^Kp,iE)dE-  jl'_^^Kp,[E)dE  (A) 

giving 


or 


AS'*-Kt'y}b  . 

Similarly, 

^A'~-  J,  KptlE)dE-  //  Kp.{E)dE  , 

£.f  ^  ^ 

(6) 


LS^^Kt^V}h 


(7) 


tion  only  becomes  pronounced  at  large  biases.  Such  a 
dependence  is  not  observed  for  filled*sute  images  because 
of  the  influence  of  the  Si  dangIing*bond  sutes. 

As  we  have  discussed  earlier,  the  empty  Si  dangling- 
bond  sutes  are  more  delocalized  or  diffused.  Their  con* 
trihution  relative  to  that  of  graphite  sutes  includes  a 
geometric  factor  of  roughly  e  where  d  is  the  work 
function  in  Ry  and  d  is  the  distance  between  the  graphite 
layer  and  the  Si  surface  layer  in  Bohr  radius.  This  factor 
has  an  order  of  l0“^  which  implies  a  negligible  contribu¬ 
tion  to  the  empty*sute  images.  On  the  other  hand,  for 
positive  lip  biases  measuring  the  filled  sutes,  the  very 
strong  localized  Si  dangling-bond  sutes  should  be  able  to 
compensate  for  the  deficiency  caused  by  the  geometric 
factor  and  contribute  significantly  to  the  filled-sute  im¬ 
ages.  This  may  be  the  reason  why  the  corrugations  m  the 
fiUed-sute  images  in  Figs.  1  and  2  are  insensitive  to 
Changes  in  tip  bias. 


V.  CONCLUSIONS 

The  ^screpancy  between  the  6X6  STM  geometry  and 
the  6v'3X6v'2  LEED  paiiem  for  the  ^-SiCdll)  surface 
can  be  explained  by  an  incommensurate  graphite  mono- 
layer  grown  on  the  Si-terminated  ^•SiCd  1  \)  surface  dur¬ 
ing  the  annealing  process.  The  contrast  reversal 
phenomenon  observed  in  the  STM  images  at  opposite  tip 
biases  can  be  explained  by  the  different  electronic  contri¬ 
butions  from  graphite  sutes  of  those  C  atoms  with  Si 
atoms  directly  (or  almost  directly)  underneath  and  those 
without.  The  Si  dangling-bond  sutes  contribute 
significantly  to  the  corrugations  observed  in  the  filled- 
sute  images. 


The  opposite  signs  in  Eqs.  (5)  and  (7)  explain  the  contrast 
reversal  observed  in  the  (a)  and  (b)  images  in  Figs.  1  and 

In  £q.  (7)  for  empty-sute  images,  the  contrast  has  a 
dependence,  which  means  the  contrast  is  large  for  Urge 
negative  tip  bias,  hut  decreases  rapidly  as  the  bias  de* 
creases.  This  may  ezpiain  the  corrugation  dependence 
shown  in  Fig.  3,  i.e.,  the  eropty*sute-  images  show  little 
eorrugation  for  small  negative  tip  bias,  and  the  corruga* 
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Relaxation  of  the  ZnTe  and  CuCI  (110)  surfaces 
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(Received  1  October  1991;  accepted  31  December  1991) 

Ab  initia  molecular  dynamics  simulations  of  the  ZnTe  and  CuQ  ( 1 10)  surfaces  are  employed  to 
study  surface  atomic  relaxation.  We  believe  that  these  are  the  first  such  computations  for 
heteropdar  semiconductors,  and  for  their  surfaces  in  particular.  The  molecular  dynamics 
follows  the  Sankey-Niklewski  method,  and  electrostatic  interactions  are  incorporated  using 
Ewald’s  scheme  for  Gaussian  atomic  charge  distributions.  Hence  the  electrons  are  treated  in  the 
local*density  approximation,  forces  are  computed  using  the  Hellmann-Feynman  method,  and 
atoms  move  to  equilibrium  according  to  Newton’s  laws.  Using  ’’dynamical  quenching,"  we 
allow  the  "ideal"  surfaces  to  relax  according  to  these  laws  of  physics  and  then  address  a 
controversy  concerning  whether  Coulomb  forces  can  play  a  significant  role  in  determining  the 
( 1 10)  zinc  blende  surface  relaxation*.  Coulomb  effects  are  not  negligible  for  ZnTe  ( 1 10)  and  are 
as  dominant  as  covalent  effects  for  CuCl  (110).  They  reduce  the  (almost  rigid)  bond  rotation 
angle  Ui  due  to  movement  of  the  anions  up  out  of  the  (110)  surface  by  about  6*  for  ZnTe,  and 
cause  the  CuCl  anion  rotation  angle  to  be  nearly  zero. 


1.  INTRODUCTION 

Low-energy  electron  diffraction  (LEED)  analysis  is  a  very 
useful  tool  for  determining  the  atomic  positions  at  a  sur- 
ftce.  This  analysis,  as  currently  practiced,'  uses  various 
methodologies  to  determine  which  atomic  geometry  yields 
electron  diffraction  intensities  that  best  fit  data.  Simplified 
models  of  solids  are  usually  used  to  predict  these  intensi¬ 
fies.  which  are  rather  difficult  to  calculate.  One  such  model 
divides  the  solid  into  ion  cores  and  a  uniformly  distributed 
dectron  gas;'  the  ion  cores  are  singly  ionized  atoms,  inde¬ 
pendent  of  the  ioniciiies  of  the  materials.  Using  this  model. 
Duke  et  at.  have  done  LEED  analyses  of  the  ZnTe  (110) 
surface^  and  concluded  that  the  relaxation  of  this  surface  is 
as  large  as  that  of  GaAs.  with  a  surface  bond-rotation 
ingle  <t>]  of  30*.  Based  on  these  and  similar  analyses  for 
other  zinc  Uende  semiconductors,  those  authors  have  hy¬ 
pothesized  that  ( 1 10)  surfaces  of  all  zinc  blende  semicon¬ 
ductors  have  virtually  the  same  surface  relaxation,  a  nearly 
rigid  rotation  of  the  anions  up  out  of  the  surface  through 
an  ianicity-independent  angle  Uj  s;  29*.  It  is  difficult  to  un¬ 
derstand  how  the  relaxation,  which  for  low-ionicity  mate- 
risk  is  dominated  by  covalent  forces,  can  be  independent 
of  ionicity,  since  the  Coulomb  forces  by  the  subsurface 
cations  should  pull  the  surface  anions  back  down  toward 
the  surface  plane;  for  large  enough  ionicity.  those  forces 
should  largely  cancel  the  relaxation  due  to  covalent  forces. 

The  present  paper  addresses  the  issue  of  whether  or  not 
such  a  universal  (110)  surface  relaxation  is  to  be  expected 
theoretically,  using  a  first-principles  local-density  approxi¬ 
mation  theory  that  neither  divides  up  the  solid  arbitrarily 
nor  replaces  the  actual  valence  charge  distribution  with 
uniform  charge  densities.  Furthermore,  the  present  theory 
allows  the  surface  to  relax  according  to  Newton’s  laws, 
with  the  interatomic  forces  calculated  quantum  mechani¬ 
cally.  No  a  priori  model  of  the  uliimaie  relaxation  geome¬ 
try  is  fed  into  the  calculations;  the  equations  of  motion,  as 
implemented  by  our  formalism,  determine  the  relaxation 
entirely. 


Most  of  the  data  analyses  which  support  the  universal 
rigid  rotation  model  have  involved  low-ionicity  semicon¬ 
ductors  such  as  GaAs.  ZnTe  is  a  notable  exception,  having 
a  relatively  large  ionicity  and  a  relatively  large  (longitudi¬ 
nal)  effective  charge.^  To  date.  LEED  analyses  of  the  CuCl 
(110)  surface  have  not  been  presented,  to  our  knowledge, 
but  this  semiconductor  has  one  of  the  largest  ionicittes  of 
all  of  the  zinc  blendes.  The  Coulomb  interactions  of  effec¬ 
tive  charges  should  tend  to  suppress  surface  relaxation, 
provided  surface  atoms  have  comparable  effective  charges 
to  those  of  bulk  atoms.  The  above-surface,  negatively 
charged,  relaxed  anions  are  attracted  back  toward  the  sur¬ 
face  by  near-surface  cations.  The  pseudofunction  theory* 
included  these  Coulomb  effects,  as  well  as  the  covalent 
physics  of  most  other  models,  and  concluded  that  ui  for 
ZnTe  (110)  should  be  about  5*  smaller  than  for  GaAs 
(110),  almost  entirely  due  to  Coulomb  interactions.  Thus 
there  is  a  controversy  between  Duke  et  at.  on  the  one  hand, 
who  propose  that  the  reconstruction  angle  <t>i  is  virtually 
independent  of  ionicity  and  Tsai  et  at.  on  the  other  who 
predict  an  ionicity  dependence.  Duke  et  at.  base  their  ar¬ 
guments  on  analyses  of  data  with  theory  that  does  not 
handle  Coulomb  effects.  Tsai  et  at.  have  purely  theoretical 
arguments.  The  theories  predict  virtually  the  same  relax¬ 
ation  for  GaAs,  differ  somewhat  for  ZnTe,  and  differ 
greatly  for  CuCl,  where  the  ionicity  is  sufficiently  large  that 
Tsai  et  at.  predict  that  the  Coulomb  forces  inhibit  relax¬ 
ation,  leading  to  a  value  of  near  zero. 

Intermediate  between  the  theories  of  Tsai  et  at.  and 
Duke  el  al.  is  pioneering  work  by  Chadi,^  who  calculated 
the  surface  atomic  geometries  using  a  total  energy  method 
based  on  (bulk)  tight-binding  parameters  and  phenomeno¬ 
logical  ideas  about  ionicity.  Chadi  concluded  that  the  ro¬ 
tation  angles  were  virtually  indepetident  of  ionicity,  al¬ 
though  the  heights  that  the  surface  anions  assume  above 
the  surface  decrease  with  ionicity. 

Neither  the  Tsai  et  al.  theory  nor  the  Duke  et  al.  work 
actually  predicts  the  positions  of  surface  atoms;  rather  both 
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^assume  i  surface  geometry.  Then  Tsai  ef  al.  predict  its 
electronic  properties,  including  total  energy,  to  verify  that 
the  assumed  geometry  is  correct,  namely  that  it  has  a  (lo¬ 
cally)  minimum  total  energy.  In  most  of  the  work  of  Duke 
ft  the  assumed  geometries  arc  varied  to  optimize 

agreement  with  selected  data,  such  as  LEED  intensities.  As 
a  rule,  only  the  surface  and  one  to  three  subsurface  layers 
are  allowed  to  relax  in  the  models  involved  in  LEED  anal¬ 
yses.  Given  this  situation,  we  thought  it  might  be  instruc¬ 
tive  to  allow  the  surface  to  relax  according  t.o  the  laws  of 
physics,  with  no  preconceived  model  of  how  the  relaxation 
occurs.  Our  approach  uses  a  “solid"  that  is  six  atomic 
layers  thick  (and  so  should  contain  all  of  the  physics  of  the 
previous  models),  and  employs  the  well-established  local- 
density  approximation  in  a  first -principles  Newtonian  mo¬ 
lecular  dynamics. 

II.  OUTLINE  OF  METHOD 

Our  approach  is  based  on  the  very-efficient,  first-princi¬ 
ples,  real-space  method  of  Sankey  and  Niklewski  for  mo¬ 
lecular  dynamics  simulations  of  covalent  materials.^  Forces 
on  atoms  are  computed  using  the  quantum-mechanical  lo¬ 
cal-density  approximation  in  a  tight-binding-like  formal¬ 
ism,  and  the  atoms  move  in  response  to  these  forces  ac¬ 
cording  to  Newton's  laws  of  motion.  This  method  has 
worked  well  for  various  systems  of  C  and  Si  atoms.^’  In 
the  original  method,  a  neutral-atom  approximation  is  used 
to  simplify  the  formalism  and  to  enhance  the  efficiency  of 
computation.  This  approximation,  however,  also  limits  the 
application  of  the  method  to  materials  with  small  effective 
charges. 

In  the  present  work,  we  have  extended  the  method  to 
more  ionic  systems  by  incorporating  Ewald’s  technique*^ 
to  calculate  the  long-ranged  Coulomb  potentials  associated 
with  non-neutral  atoms.  This  scheme,  which  is  presented 
in  detail  elsewhere,'*  is  computationally  very  fast,  and  so 
we  arc  able  to  extend  the  multiccnter  light-binding  molec¬ 
ular  dynamics  method  to  a  wide  range  of  materials,  with¬ 
out  trading  off  much  computational  efficiency. 

The  essence  of  the  new  treatment  of  charge  transfer  is 
that  a  Gaussian  charge  distribution  centered  on  each 
atomic  site,  with  its  charge  chosen  to  neutralize  the  atomic 
charge,  is  both  added  and  subtracted.  Then  the  atoms  with 
their  neutralizing  Gaussian  charges  are  treated  as  the  neu¬ 
tral  pseudoatoms  of  the  Sankey-Niklewski  method,  and 
the  negative  of  the  array  of  neutralizing  Gaussian  charges 
is  treated  using  Ewald's  method.  (Note  that  overall  the 
array  of  Gaussian  charges  is  neutral.)  The  numbers  of  s 
and  p  valence  electrons,  as  defined  by  the  occupation  num¬ 
bers  of  5  and  p  orbitals  for  a  given  atom,  are  determined 
self-con sistenily,  and,  in  general,  arc  nonintegers.  The  de¬ 
tailed  algebra  of  these  modifications  of  the  neutral-atom 
theory  are  given  in  Ref  11. 

The  ZnTc  (110)  surface  is  modeled  by  a  six  layer  slab 
repeated-slab  model  (i.c.,  a  super-cell  model)  with  periodic 
boundary  conditions  in  three  dimensions.  The  vacuum  re¬ 
gion  between  adjacent  slabs  is  8.5  A  thick.  The  ideal  sur¬ 
face  plane  corresponds  to  r  »  0,  and  the  relaxed  anions  lie 
in  the  upper  half-plane  for  ZnTe  {<o\  >  0).  The  Observed 
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lattice  constant  of  bulk  ZnTe,  6.101  is  used  to  deter¬ 
mine  the  reference  lattice  and  to  define  basts  and  redpro- 
caMattice  vectors.  (We  have  tested  that  an  adjustment  of 
the  lattice  consunt  by  0, 1  k  does  not  significantly  alter  the 
surface-bond  roution.)  For  CuCI,  we  treat  the  Cu  d  elec¬ 
trons  as  core  electrons  rather  than  valence  electrons,  and 
this  causes  the  theoretical  bulk  lattice  constant  to  differ 
somewhat  from  the  experimental  one.  Hence,  we  use  the 
theoretical  value,  4.70  A  here,  in  order  to  be  compleldv 
consistent  theoretically.'-'  The  masses  of  atoms  in  the  bo(- 
tom  surface  of  the  slab  are  set  to  huge  numbers,  so  thit.  for 
all  practical  purposes,  these  atoms  are  motionless  and  arc 
regarded  as  fixed  bulk  atoms.  All  atoms  in  other  layen 
have  atomic  masses  (66.38,  127.60.  63.546,  and  35.45)) 
a.u.  for  Zn,  Te,  Cu,  and  Cl,  respectively.  The  unit  cell 
contains  twelve  atoms,  c.g.,  one  Zn  and  one  Tc  atom  pci 
layer.  We  use  s  and  p  orbitals  centered  on  all  aioms  ind 
treat  the  d  orbital  of  Cu  as  a  core  electron  that  is  contained 
in  the  pseudopotential  and  does  not  participate  in  the  mo¬ 
lecular  dynamics.  Molecular  dynamics  simulations  arc 
started  assuming  the  ideal,  unrelaxed  zinc  blende  (110) 
surface.  Movements  of  all  atoms  arc  governed  by  the  total 
forces  acting  on  them.  During  the  simulation,  the  velocities 
of  all  atoms  are  quenched  (set  to  zero)  when  the  insun 
taneous  temperature  for  a  given  time  step  is  smaller  than 
for  the  previous  time  step.  This  means  that,  overall,  the 
atoms  decelerate.  Each  time  step  is  2.074  fs  long.  The  sim¬ 
ulation  is  terminated  when  all  atomic  positions  subilizz 
For  each  time  step,  the  charges  or  the  occupation  nurhben 
of  the  valence  s  and  p  orbitals  on  all  atoms  are  determined 
self-consistently,  within  a  tolerance  of  0.01  electrons.  The 
self-consistency  in  charge  is  more  readily  achieved 
each  successive  time-step. 

111.  RESULTS 

Using  the  original  zero-charge-transf  er  method  of  San- 
key  and  Niklewski,^  we  obuin  oi  ^  iV  and  a  second-layer 
bond- rotation  angle  cuj  5*  for  ZnTc  (110).  The  corre 
spending  results  for  CuQ  arc  »  +31*  and 
3=  —  T.  Tlic  results  for  ZnTc  and  CuCI  arc  similar  to  the 
well-known  surface  relaxation  ofGaAsfnO):  »  30.f 
and  0)2  =  —  6.2*  by  this  method.  GaAs  has  a  small  effec¬ 
tive  charge,  and  so  is  close  to  the  purely  covalent  limii. 
stated  previously,  the  Sankcy-Niklcwski  method  assumes 
that  all  atoms  are  neutral  and  there  is  no  long-ranged  Cou¬ 
lomb  potential.  Thus  these  zero-charge-transfer  results  an 
be  regarded  as  the  covalent  limit  in  which  the  surface  an- 
ion.s  rise  up  out  of  the  surface  and  the  cation.s  descend 
toward  the  bulk,  such  that  the  surface-bond  rotates  b\ 
about  29*±2*  (the  commonly  accepted  value  of  0)^). 

With  the  inclusion  of  charge- transfer  and  long-ranged 
Coulomb  forces,  we  find  that  the  relaxation  of  the  ZnTr 
(110)  surfaces  still  can  be  described  as  a  surface  bond 
rotation  and  a  second  layer  counter  rotation.  (For  CuCI. a 
similar  model,  but  w*ith  more  bond  contraction  and  distor¬ 
tion.  is  necessary.)  Hence  the  model  is  described  by  ihf 
two  angles  a)|..and  0)2.  and  the  2  displacements  of  anions 
and  cations.  (Sec  Figs.  1  and  2.)  The  displacemenis  nfihe 
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Fig.  1.  Tbc  vuw  of  the  reiexed  zinc  blende  (110)  surface.  Smaller 
sad  larger  circles  correspond  lo  cations  and  anions,  respectively.  The 
ingle  between  ihe  plane  (the  unrelaxed  surface)  and  the  surfaca 
lAKMis  is  a>|.  The  dark  circles  are  esseniially  infinitely  massive  in  the 
preseni  model. 


near*$urface  atoms  predicted  by  the  Simula)  ions  are  given 
in  Tables  I  and  II  for  ZnTe  and  CuCl,  respectively. 

In  the  case  of  ZnTe,  )he  predicied  bond-rotation  angle 
u)i  with  charge  transfer  is  signihcanily  reduced  in  compar¬ 
ison  with  both  GaAs  and  zero-charge-transfer  ZnTe,  being 
rjj  »  24.6*;  in  this  case  we  have  ^  4.6*.  Hence  o>|  is  about 
6*  smaller  than  obtained  for  neutral  aioms  with  zero- 
charge  transfer.  (This  compares  favorably  wuh  a  corre¬ 
sponding  difference  of  5*  obtained  using  the  pseudofunc- 
lion  method,  fixing  the  bond  lengths,  and  allowing  only  ihe 
top  layer  of  atoms  to  relax. 

For  CuCl,  there  is  almost  no  surface  bond  rotation;  )he 
theory  predicts  oii  =  -i-  3.6*  and  ^  6.2*.  The  mo¬ 

lecular  dynamics  calculations  that  we  have  executed  as¬ 
sume  that  the  Cu  d  electrons  are  part  of  the  core.  Hence 
Ihe  coupling  of  the  Cu  d  electrons  to  the  sp  valence  band 
has  been  neglecied. 

In  order  to  determine  that  this  neglect  of  dsp  coupling 
is  a  reasonable  approximation,  we  have  done  pseudofunc- 
lion  calculations,^^  which  do  account  for  the  coupling,  for 
a  four-layer  isolated  slab.  The  pseudofunction  calculations 
are  compuiaiionally  intensive,  so  the  only  relaxaiton  that  is 
allowed  is  a  rigid  rotation  of  the  surface  layer  through  an 
angle  When  we  do  these  calculations  with  the  d  elec- 


FiC.  2.  lllusiraiex  the  coordinate  syitem  with  x  positive  in  the  |00l] 
duictiOA,  y  poaiiivc  in  ihe  (llO]  direction,  and  x  positive  normal  to  the 
wneiaxad  ( 1 10 )  surfaca. 


Table  I.  Calculated  atomic  dtaplacement  M  (in  A)  with  rmpect  to  the 
idea)  truncated  bulk  fcometry,  as  calculated  for  the  ZnTe  ( tlO)  surface; 
dit  m  (dx,4r,dz)  where  x  ia^ normal  to  the  aurfacc.  a  b  in  the  [(X31] 
direction,  and  is  in  the  {iTo]  direction.  Coordinates  for  fourth  and 
deeper  layen  should  be  viewed  as  unrehable,  because  the  model  fixes  the 
bottom  (sixth)  layer. 


Layer,  atom  AR 

Pint,  2ii(  0  27.  -  0.04,  -  0.27) 
First,  Te(  h-  0. 13.  -  0.02,  +  0  37) 
Second.  Zn(  -  0.38,  +  0.01,  H-  0. 19) 
Second,  Te(  -  0.03,  0.04.  0.10) 

Third.  2n(  +  000.  +  0.00,  +  0.20) 
Third.  Te(  +  0Xi&,  +  0  00.  -  0.01 ) 

Founh.  2n(  -O.Ol,  +  0  00.  -  0.21) 
Founh,  Te(  +003,  -r  0.00, 0.13) 
Fifth.  2n(  -0.24.  +  0.00,  -t-  0.02) 
Fifth.  Te(  -  0  1 1 ,  -r  0.00.  -  0  01) 
Sixth.  Zn(  +  0.00,  -r  0  00.  +  0.00) 
Sixth.  Te(  +  0  00.  0.00,  ^  0.00) 


irons  consirained  to  the  core,  we  find  ^  —0.1 3*.  only 
3.T  different  from  the  comparable  molecular  dynamics  re* 
suit.  When  we  treat  the  Cu  d  electrons  as  valence  electrons 
and  let  them  couple  to  O’s  s  and  p  electrons,  we  find  <u, 
=s  -(■  1.44*.  comparable  with  both  the  molecular  dynamics 
and  the  pseudofunciion  results  for  d  electrons  in  the  Cu 
core.  Thus  we  conclude  that  our  treatment  of  the  Cu  d 
electrons  introduces  an  insignificant  uncertainty  to  the  re¬ 
sult  for  (i)|  of  about  1.6*. 

Tables  I  and  II  give  our  predicted  relaxed-surface  coor¬ 
dinates  for  ZnTe  (110)  and  CuCl  (110).  The  final  config¬ 
urations  were  obtained  by  the  surface  atoms  successively 
distorting  and  accelerating  from  the  “ideal”  surface  m  time 
intervals  of  2.074  fs.  After  a  time  interval  in  which  the 
kinetic  energy  of  the  surface  increases,  the  surface  is 
“cooled”  by  having  all  its  atomic  velocities  set  to  zero. 
Wltile  this  is  a  kinder  and  gentler  way  of  relaxing  to  equi¬ 
librium  than  cleavage,  it  nevertheless  should  provide*  a 
good  model  of  the  (110)  surface's  equilibrium  state.  There¬ 
fore,  we  believe  that  the  qualitative  conclusion  that  <U|  de- 

TaBLC  It.  Calculaied  aiotnic  dUptMcment  (in  k.)  wiih  ictpeci  lo  ihe 
ideal  truncated  bulk  geomeiry,  as  calculated  for  the  CuCl  ( IIO)  surface: 
4R  —  where  x  is  normal  to  the.  surface,  x  it  in  the  |00I) 

direction,  end  y  is  in  the  (iTO]  direction. 


Layer,  atom  aR 


First,  Cu{  +  0.32,  -  0. 12,  -  0. :  I ) 
First,  a(  -021.  +  017. -0.17) 
Second.  Cu(  +0.16,  -^0. 10. +  0.13) 
Second.  CI(  -0.34.  -  0.20  -  0.05) 
Third,  Cu(  +0.21. -on. +0.10) 
Third.  a(  -0.28, +  0.19, -0.13) 

Fourth.  Cu(  +  0.24. -r  0.13.-^005) 
Founh.  CI(  -0.30,-018.-0  07) 
Fifth.  Cu<  +O.I9.  +  0  13.-0.I4) 
Fifth.  Cl(  -026,-0.20.-015) 
Sixth.  Cut  +  0  00.  +  0.00.  T  0  00) 
Sixth.  Cl(  •!-  0  00. 1-  0  00.  +  0.00 ) 
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cratsts  sfgnificantly  with  ionicity  must  be  correct  We  be¬ 
lieve  ihat  only  if  the  k)Cal-den»ty  approximation  for 
computing  atomic  forces  is  inadequate  for  computing 
charge  transfer/’  or  if  our  slab  is  too  thin,  could  ihese 
conclusions  ble  reversed.  Since  the  local-density  approxima¬ 
tion  has  a  long  history  of  being  adequate,  and  bemuse  we 
have  studied  the  dependence  on  super-cell  size  (See  Ap¬ 
pendix),  we  think  that  future  progress  on  this  problem  is 
likely  to  come  from  reanalysis  of  LEED  data  and  from 
LEED  experiments  on  CuCl  (110). 

From  our  results,  the  electrostatic  effects  on  zinc  blende 
(110)  surface  relaxation  are  clearly  demonstrated.  Our  re¬ 
sults  agree  with  previous  calculations  of  Tsai  e/  al.  using 
the  firsi-principles  pseudofunction  method,  and  disagree 
with  ihe  analyses  of  LEED  intensities. 

IV,  SUMMARY 

In  summary,  using  an  improved  ab  iniiio  multicenter 
reaLspace  tight*binding*like  method  of  molecular  dynam¬ 
ics  simulations  that  allow*s  charge  transfer  among  atoms, 
we  found  thai  the  relaxations  of  the  ZnTe  ( 1 10)  and  CuCl 
(110)  surfaces  are  significantly  smaller  than  that  of  GaAs 
(110),  due  to  Coulomb  interactions  between  atoms.  In  the 
case  of  CuCl  the  effect  was  so  dramatic  that  there  was 
almost  no  relaxation  at  all.  Our  results  imply  that  a  com* 
plete  understanding  of  the  zinc  blende  (110)  surfaces  re¬ 
quires  an  adequate  treatmeni  of  charge  transfer. 
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APPENDIX 

First  principles  molecular  dynamics  calculations  of  the 
sort  presented  here  are  computationally  intensive,  and 
hence  demand  a  limited  size  for  the  super  cell.  We  have 
used  a  rather  large  slab  with  six  layers,  essentially  fixing 
the  positions  of  the  aloms  in  the  bottom  (sixth)  layer.  If 
the  atomic  displacements  from  lattice  sites  in  the  fiflh  layer 
were  nearly  zero,  we  could  be  completely  confident  that 
our  results  had  fully  converged  with  respect  to  slab  size. 
However  ihe  atomic  displacements  in  the  ihird  to  fifth 
layers  are  significant  (see  Table  I) — a  reflection  of  the  fact 
that  elastic  forces  are  long  ranged,  and  a  source  of  concern 
thai  the  slab  may  be  too  thin. 


Ta8U  111.  Cafettiatcd  m  Mbomk  diiplarmfafi  Alt  im  A) 

ZiiTe  ( 1 10)  fyrface  Ibr  a  iU4tycr  aiodd  nirnm  a  Avc4aycr  moM;  gAK 
V  g(Ax.Ay.Ai)  where  X  is  aomal  to  the  aaitee,  jr  is  m  the  (001]  4sf«e* 
(ton  and  y  is  in  the  (iTO]  dtrectiofi. 

Layer,  alom  gAR 

Fim.  Zn(  ^  0  02. 4*  0.05,  +  0.05) 

Fim.  Te(  4-  0.07.  4*  0.05.  +  0.03) 

Second.  Zn(  -  003,  ^  0.06.  -  002) 

Second.  Te(  +  001. 4*  0  05.  +  0.06) 


To  address  the  dependence  on  slab  size  within  the  strin 
gent  limiutions  imposed  on  us  by  available  computational 
resources,  we  computed  a  five-layer  model  of  ZnTe  (tlO) 
and  compared  it  with  the  corresponding  six-layer  model. 
The  differences  in  the  angles  and  toj  were  ^  1.4*  and 
—  3*,  respectively.  The  corresponding  differences  in  dis 
placements  are  given  in  Table  111,  and  are  a  good  indict 
tion  of  uncertainties  in  the  calculations. 
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'^With  Ihe  cxcepifon  of  Table  111.  the  parameter  of  ihe  Sankn* 
Niklewski  model  is  2.38  A  for  Cu,  3.24  A  for  Te,  2.78  A  for  Zn.  sad 
2.65  A  for  Cl.  For  Table  Ill.  we  used  %  »  2.12  A  for  Zn  and  3.06  k  far 
Te.  The  parameter  r,  is  ch(mn  moally  for  convenience,  provided  i 
reasoiiabk  charge  Iranffer  can  be  obcainad.  as  determined  b>  f/% 
where  ^  is  i he  1on|Uudinal  dTactivc  charge  and  ^  k  Ihe  static  didacinc 
consuni.  The  bulk  lattice  consUni  of  CttO  is  5.4057  A. 

'^R.  V.  Ktsomaki.  M*H.  Tsai.  T.  N.  Rhodin,  and  D.  D.  OiamMiss.  PKn 
Re%.  B34.  2656  (1986). 

'^The  charge  iran^fcr  obtained  here  is  eaaenllally  the  same  as  obtained  m 
Ref  4  We  did  note  some  sensitivity  of  ihc  calculalad  charge  irtnsfcrie 
Ihc  cutoff  radius  v  Hence  we  beltcvt  that  the  charge  iransTer  of  Rc(  4 
is  (he  charaeicfistic  value  10  be  expected  from  localnleiitiiy  ihcory. 
independeni  of  how  that  theory  is  impkincniad  compuiaiiona%. 
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Using  a  scanning  tunneling  microscope,  quantum  dots  with  diameters  of 
order  25  A  have  been  fabricated  on  the  cleaved  InSb(llO)  surface  in 
UHV.  Both  In  and  Sb  atoms  were  clearly  resolved  on  the  cleaved  surface. 


Scanning  tunneling  microscopes  (STM's)  present 
the  opportunity  not  only  to  image  a 
semiconductor  surface  with  unprecedented 
resolution  but  also  to  controllably  alter  the 
surface  on  the  Angstrom  scale,  as  a  result  of 
the  tip* surface  interaction,  Becker  et  al, 
demonstrated  that  single  atoms  could  be 
manipulated  on  the  Ge(lll)  surface  with  a  single 
voltage  pulse  [1],  A  number  of  authors  [2]  have 
"machined"  various  surfaces'  on  the  nanometer 
scale,  either  by  crashing  the  STM  tip  into  the 
surface,  by  applying  voltage  pulses,  or  by  local 
heating  **  with  experiments  performed  in  air, 
under  liquids,  or  in  ultra*high  vacuum.  In  this 
paper  we  present  some  of  the  first  [3]  STM 
images  of  the  InSb(llO)  surface  under  ultra* high 
vacuum  and  we  show  how  quantum  dots  can  be 
controllably  engraved  in  the  surface. 

Ve  employed  a  Pachyderm*4  sculpted  STM  (^1( 
which  is  machined  almost  in  its  entirety  from  a 
single  block  of  stainless  steel  and  hence  has 
unusual  \ibratlonal  stability.  With  this  STM,  we 
obtained  the  Images  of  Fig.  1  and  Fig.  2  by 
scanning  the  InSb(llO)  surface.  This  surface  was 
cleaved  in  10*^^  Torr  vacuum,,  and  the  images 
were  taken  with  a  tunneling  current  of  *120  pA. 
Fig.  1  was  taken  with  a  negative  sample  bias  of 
*0.2  V,  and  so  is  sensitive  to  the 
electronically  occupied,  Sb*like  states  of  the 
valence  band.  The  prominent  features  of  the 
image  are  Sb  atoms,  spaced  by  the  accepted 
lattice  constants  (assuming  a  29”  Rigid  Rotation 
Model  15]  of  the  surface  relaxation).  The 
number  of  visible  defects  in  this  region  of  the 
surface  is  about  one  per  *10^  atoms. 

Fig,  2  shows  InSb  before  and  after  a  quantum- 
dot  hole  is  fabricated  on  its  (110)  surface. 
Positive  sample  bias  of  +0.3  V  produced  the 


images  of  Fig.  2a  and  2b,  which  are  sensitive  to 
electronically  empty.  In* derived  states. 
Negative  bias  of  -0.3  V,  yielded  Fig.  2c. 

By  positioning  the  STM  tip  over  a  spot  for  * 
lO^s  we  formed  quantum  dots,  such  as  those  of 
Figs,  2b  and  2c.  The  positions  of  the  dots  were 
easily  controllable.  The  perfect  InSb(llO) 


Fig.  1,  STM  image  of  InSb(llO)  with  atomic 
resolution.  The  surface  was  cleaved  in 
ultra-high  vacuum  (10*  torr  range)  and 
imaged  at  *0.2  V  sample  bias  with  120  pA 
tunneling  current.  With  negative  bias  the 
Sb  atoms  are  prominent  and  the  In  atoms  are 
slmost  invisible.  The  data  were  Fourier 
filtered. 
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surface  was  more  susceptible  to  disruption  by 
the  STM  tip  than  other  surfaces  we  have  studied 
[2,6),  and  so  preservation  of  the  perfect 
surface  required  rapid  scanning,  while  the 
formation  of  quantum  dots  was  almost  effortless. 
The  diameter  of  the  doc  in  Fig.  2  is  about  “  25 
A  as  can  be  deduced  by  recognizing  chat  the  Sb 
atoms  of  Fig.  2b  exhibit  a  rectangular  surface 
unit  cell  of  6.5  A  X  4.6  A.  Close  Inspection  of 
Figs.  2b  and  2c  reveals  images  of  Individual  Sb 
and  In  atoms. 

The  quantum  dot  is  a  small  surface  void,  with 
a  depth  of  one  or  two  atoms,  2  A  to  4  A.  as 
determined  by  scanning  in  the  dot  The  size  of 


the  voids  can  be  increased  by  holding  the  rip 
fixed  for  a  longer  time  or  by  scanning  over  a 
very  small  area. 

These  results  demonstrate  that  Angstrom* scale 
ST.M  lithography  of  InSb(llO)  surfaces  will  be 
feasible,  and  that  it  might  be  possible,  if  Che 
electronic  structures  of  the  quantum  dots  are 
favorable,  to  fabricate  Angstrom- scale  niemory 
bits  on  this  surface. 
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Fig.  2.  STM  image  of  InSb(llO)  with  a 
sample  bias  of  (a)  +0.3  V,  (b)  +0.3  V,  and 
(c)  -0.3  V;  the  tunnel  current  was  -120  pA. 
Individual  In  atoms  are  clearly  resolved  ir. 
(a)  and  (b) ;  Sb  atoms  are  resolved  in  (cj 
A  quantum  dot  is  visible  ir.  (b)  and  (c). 
The  data  were  Fourier  filtered. 
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Using  scanning  tunneling  microscopy  (STM),  we  have  observed  an  antiphase 
boundary  that  resembles  a  perfect  dislocation  on  a  Si(111)-(7x7)  surface.  The  extra 
row  of  adatoms  forming  the  boundary  lies  on  the  untauited  hafi  of  the  7x7  unit  cell,  in 
agreement  with  total  energy  calculations  using  the  first-principles  self-consistent 
pseudofunction  method.  The  filled-state  STM  Image  shows  a  missing  Interior  adatom 
on  the  unfaulted  half,  in  agreement  with  calculations  of  the  partial  density  of  states  of 
an  adatom  surrounded  by  three  rest-atoms. 


PADS  numbers:  68.35.Bs;  61.16Di;  61.70.Jc. 


In  recent  years,  scanning  tunneling  microscopy  (STM)  has  been  applied 
extensively  to  the  study  of  surface  defects.  Specifically,  atom-resolved  images  of 
dislocations  have  been  obtained  on  surfaces  of  PbS(00l)\  Au(lll)^,  GaAs(llO)®, 
and  PtNi(ll1)^.  On  the  well-studied  Si(lll)-(7x7)  surface,  however,  no  dislocations 
have  ever  been  found,  although  domain  boundahes  are  frequently  observed^'^. 
These  boundahes  generally  separate  two  adjacent  7x7  domains  on  the  same  terrace, 
with  a  phase  slip  across  the  boundary.  The  slippage  is  usually  taken  up  by  missing 
adatoms,  and  the  boundary,  more  often  than  not,  follows  a  wavy  path.  We  present 
here  the  first  STM  images  of  an  antiphase  boundary  that  resembles  a  perfect 
dislocation  on  a  Si(111)-(7x7)  surface.  We  will  also  compare  the  filled-state  image  of 
this  boundary  with  first  principles  calculations  of  the  electronic  structure  in  the  vidnity 
of  the  boundary. 

The  sample  used.in  the  experiment  was  a  n-type  Si(l  1 1 )  wafer  with  a  resistivity 
of  0.015  G  cm.  The  surface  was  cleaned  the  conventional  way  by  flashing  it  to 
1200^0  for  several  seconds.  The  base  pressure  of  the  ultrahigh  vacuum  chamber  was 
51x10*’°  torr  and  rose  to  2x10*°  torr  during  flashing.  After  cooling  to  room 
temperature,  STM  images  showed  the  familiar  7x7  reconstruction.  With  the  sample 
bias  at  2V,  Le.  sampling  unoccupied  states,  we  observe  an  extra  row  of  adatoms 
bisecting  the  two  halves  of  a  7x7  unit  cell  as  indicated  by  arrows  in  Fig.  1(a).  The 
filled-state  image  of  the  same  area  is  shown  in  Fig.  1(b),  and  we  observe  that  one  of 
the  middle  adatoms  is  missing  on  the  unfautted  half  of  the  7x7  unit  cell.  The  extra  row 
of  adatoms  spans  the  entire  width,  ~350A,  of  a  7x7  terrace,  running  from  an  upstep 
edge  as  shown  in  Fig.  2(a),  to  a  downstep  edge  in  Rg.  2(b). 

Because  of  Its  uninterrupted  straight-line  nature,  this  row  of  adatoms  resembles 
a  perfect  dislocation  emerging  from  the  bulk  to  the  surface.  On  closer  inspection, 
however,  it  becomes  clear  that  the  formation  of  this  defect  does  not  require  an  extra 
plane  of  atoms  in  the  bulk  as  a  dislocation  would.  Fig.  3  depicts  a  plan  view  of  the  unit 
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cell  containing  the  defect  and  also  a  side  view  across  a  section  labelled  AB.-  From  the 
plan  view,  we  see  that  the  unfaulted  half  of  the  7x7  unit  cell  undergoes  a  phase  shift  by 
half  the  distance  between  two  adatoms,  or  3.84A,  across  the  line  defect.  However,  the 
extra  row  defect  itself  forms  the  unfautted  half  of  a  7x7  unit  cell  enclosed  by  dashed 
lines  in  the  plan  view  of  Fig.  3.  Moreover,  the  entire  unfaulted  half  (including  the  extra 
row)  almost  forms  half  of  a  9x9  unit  cell  (minus  a  comer  adatom).  From  the  side  view 
of  Fig.  3,  we  see  that  the  existence  of  an  extra  row  of  adatoms  on  the  surface  does  not 
disturb  the  arrangement  of  bulk  atoms.  So  this  particular  line  defect  can  be 
considered  simply  as  an  antiphase  boundary  on  the  7x7  surface. 

The  two  most  interesting  features  of  our  STM  observations  of  the  antiphase 
boundary  are:  (a)  the  extra  row  of  adatoms  forms  an  extension  of  the  unfaulted  half  of 
the  7x7,  i.e.  a  quasi*9x9  half  unit-cell;  and  (b)  the  middle  adatom  of  the  unfaulted  7x7, 
or  the  Interior  adatom  of  the  unfaulted  9x9,  is  missing  in  the  filled-state  image.  In  order 
to  provide  an  explanation  for  our  observations,  we  perform  total  energy,  or  zero- 
temperature  free  energy,  calculations  using  the  first-principles  self-consistent 
pseudofunction  method^. 

In  the  model  for  the  antiphase  boundary  in  Fig.  3.  the  only  adatom  surrounded 
by  three  first-layer  rest-atoms  (located  at  the  apexes  of  the  dotted  triangle)  is  the 
interior  adatom  on  the  quasi-9x9  unfaulted  half.  All  other  adatoms  in  Fig.  3,  whether 
on  faulted  or  unfaulted  half,  have  either  one  or  two  nearby  rest-atoms.  These  rest- 
atoms  are  not  directly  bonded  to  any  adatom.  In  view  of  the  demanding  task  to 
perform  calculations  on  the  7x7  surface  as  demonstrated  recently®*®,  we  simulate  the 
antiphase  boundary  (or  more  precisely,  the  quasi-9x9  half  unit-cell)  by  three-fold 
rotationally  symmetric  two-dimensional  unit  cells  with  the  adatom  at  the  origin  as 
shown  in  Fig.  4.  In  this  structure,  each  adatom  is  surrounded  by  three  rest-atoms  as 
outlined  by  dotted  lines  in  Fig.  4(a}.  Each  unit  cell  contains  one  Si  adatom  and  four  Si 
atoms  in  each  layer  of  a  six-layer  isolated  slab  model.  Four  hydrogen  atoms  are 
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attached  to  the  bottom  layer  to  saturate  the  dangling  bonds.  The  adatoms  form  a  2x2 
hexagonal  geometry.  The  difference  between  unfaulted  and  faulted  surfaces  lies  in 

the  positions  of  the  first  layer  Si  atoms  as  shown  in  Figs.  4(a)  and  4(b)  respectively. 
The  adatom  is  located  at  the  T^  site,  which  is  directly  above  a  second-layer  atom.  The 

height  of  the  adatom  above  the  second-layer  atom  is  taken  to  be  equal  to  the  bulk  SI 
bond  length,  i.e.  2.35A.  The  heights  of  the  adatoms  on  both  unfaulted  and  faulted 
halves  are  taken  to  be  the  same.  These  are  reasonable  approximations  because  a 
change  in  the  Si-Si  bond  length  by  O.OSA  only  changes  the  bond  energy  by  about 
0.03  eV^°. 

In  the  total-energy  calculations,  we  use  three  two-dimensional  special  k-points 
of  Cunningham”'^  for  a  hexagonal  lattice  to  approximate  the  integration  over  the  first 
Brillonin  zone.  Our  calculations  show  that  the  unfaulted  surface  of  Fig.  4(a)  is  more 
stable  than  the  faulted  surface  of  Fig.  4(b)  by  0.33  eV  per  atom.  The  difference  in 
energy  accounts  for  the  occurrence  of  the  antiphase  boundary  on  the  unfaulted  half, 
extending  the  7x7  to  a  quasi-9x9  geometry  on  that  half. 

We  have  also  calculated  the  partial  density  of  states  of  the  adatoms  for  the 
unfaulted  surface  shown  in  Fig.  4(a).  The  missing  interior  adatom  in  our  filled-state 
image  of  Fig.  1(b)  is  similar  to  that  reported  earlier  by  Becker  et  al.^^  for  a  filled-state 
image  of  a  Si(1 1 1)-(9x9)  reconstruction  in  v.'hich  the  interior  adatom  was  considerably 
darker  than  the  corner  and  middle  adatoms.  Becker  et  al.^^  speculated  that  because 
three  rest-atoms  surrounded  the  interior  adatom,  more  charges  were  transferred  from 
this  adatom  to  the  rest-atoms  than  either  the  corner  or  middle  adatoms  which  had 
fewer  nearby  rest-atoms.  This  increased  charge-transfer  rendered  the  Interior  adatom 
darker  than  the  rest  in  the  filled-state  image. 

Our  calculations  on  the  surface  shown  in  Fig.  4(a),  based  on  sampling  energy 
levels  of  six  special  k-points  of  Cunningham^  \  are  plotted  in  Fig.  5.  The  energies  are 
plotted  relative  to  the  calculated  Fermi  level  of  this  system,  which  is  5.03  eV  below  the 
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vacuum  level.  The  partial-den sity-of-states  curve  shows  four  peaks  located  at  -2.7  eV, 
•1.2  eV,  -0.2  eV  and  2.0  eV.  The  (dI/dV)/(W)-versus-energy  curve  for  corner  adatoms 
measured  by  Wolkow  and  Avouris'^  is  also  reproduced  in  Fig.  5  for  comparison.  The 
energy  peaks  of  the  (dI/dV)/(IA/)  curve  occur  at  -1.9  eV,  -0.4  eV,  0.5  eV  and  1.6  eV. 
Our  calculated  peaks  differ  from  the  measured  peaks  by  -0.8  eV,  -0.8  eV,  -0.7  eV  and 
0.4  eV  respectively.  The  shift  of  density  of  states  towards  lower  energies  in  the  filled 
band  corresponds  to  the  addition  of  a  positive  potential  to  the  adatom.  This  fits  the 
description  of  the  adatom  having  more  electronic  charges  transferred  to  the  three  rest- 
atoms  such  that  it  becomes  more  electropositive.  The  charge-transfer  contributes  to 
the  dark  interior  adatom  in  the  filled-state  STM  image  in  Fig.  1(b).  The  empty-state 
image  in  Fig.  1(a),  on  the  other  hand,  shows  all  equally  bright  adatoms  despite  the  fact 
that  our  calculations  show  an  upward  shift  in  energies  by  0.4  eV  in  the  conduction 
band.  This  is  probably  due  to  the  more  delocalized  or  diffused  nature  of  the  empty 
states  further  away  from  the  Fermi  level  and  the  sampling  of  these  states  is  not  affected 
by  small  shifts  in  the  partial  density. 

In  conclusion,  we  have  observed  an  antiphase  boundary  that  resembles  a 
perfect  dislocation  on  the  Si(111)-(7x7)  surface.  By  performing  calculations  on  the 
total  energies  and  partial  density  of  states  using  the  pseudofunction  method,  we  have 
explained  the  preferred  location  of  the  antiphase  boundary  on  the  unfautted  half  of  the 
7x7  and  the  missing  interior  adatom  on  the  resulting  quasi-9x9  structure  in  the  filled- 
state  STM  image. 

This  work  was  supported  by  the  U.S.  Army  Research  Office  contracts  DAAL  03- 
92-G-0038  and  DAAL  03-91 -G-0054. 
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FIGURE  CAPTIONS 

Fig,  1  (a)  Empty-state  and  (b)  fiiled-state  STM  images  of  a  disiocation-like  antiphase 
boundary,  indicated  by  arrows,  on  a  Si(11 1)-(7x7)  surface.  Scan  area  =  150A  x 
150A.  Tunnel  current  »  1  nA.  Sample  bias  =  2V  for  (a)  and  -IV  for  (b).  The  7x7 
unit  cell  is  outlined. 

Fig.  2  STM  images  of  the  antiphase  boundary  (indicated  by  arrow)  running  across  the 
full  width.  ~350A,  of  a  7x7  terrace  from  (a)  the  upstep  edge  to  (b)  the  downstep 
edge.  Scan  area  «  200A  x  200A.  Tunnel  current «  1  nA.  Sample  bias  =  2V. 

Fig.  3  A  structural  model  for  the  antiphase  boundary  showing  a  plan  view  and  a  side 
view  across  a  section  AB.  The  faulted  and  unfaulted  halves  are  labeled  F  and 
U  respectively.  The  dashed  lines  enclose  the  antiphase  boundary  forming  an 
almost  perfect  7x7  geometry.  The  unfaulted  half  has  a  quasl-9x9  structure 
(minus  a  corner  adatom)  with  an  interior  adatom  surrounded  by  three  rest- 
atoms  at  the  apexes  of  the  dotted  triangle.  This  interior  adatom  appears 
missing  in  the  fiiled-state  image  in  Rg.  1(b). 

Fig.  4  Plan-view  models  of  the  unfaulted  surface  (a)  and  faulted  surface  (b)  used  in 
our  calculations.  •  Si  adatom;  O  first-layer  Si;  o  second-layer  and  third-layer 
Si;  x  fourth-layer  Si.  in  this  model,  each  adatom  is  surrounded  by  three  rest- 
atoms  located  at  the  apexes  of  the  dotted  triangle  in  (a).  The  Si  adatoms  form  a 
2x2  unit  cell  outlined  in  both  (a)  and  (b). 

Fig.  5  (a)  The  partial  density  of  adatom  states  calculated  for  the  surface  shown  in  Fig. 
4(a).  (b)  The  (dI/dV)/(17V)  curve  of  Wolkow  and  Avouris''^  for  a  corner  adatom 
on  a  7x7  surface.  The  energy  peaks  in  each  curve  are  labeled  1  through  4. 
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Scanning  tunneling  microscopy  (STM)  images  of  the  (110)  surface  of  InSb  reveal 
apparent  surface  interstitial  Sb  atoms  when  the  sample  bias  voltage  is  small  and 
negative,  corresponding  to  tunneling  of  electrons  from  the  top  of  the  valence  band 
of  InSb.  These  surface  “interstitials”  disappear  at  other  bias  voltages.  The  ob¬ 
servations  of  these  “ethereal  interstitials”  are  explained  in  terms  of  below- surface, 
second-layer  Sb  atoms,  that  appear  to  lie  above  the  surface  due  to  the  re-hybridized 
character  of  the  surface  Sb  wavefunctions  near  the  top  of  the  valence  band. 


1.  Introduction 


In  this  Communication,  we  report  scanning  tunnel¬ 
ing  microscopy  (STM)  measurements  of  the  InSb  (110) 
surface  under  conditions  of  small  negative  sample  bias 
(tunneling  of  electrons  from  near  the  top  of  the  InSb 
valence  band  to  the  STM  tip).  We  find  an  apparent 
face-centered  rectangular  surface  crystal  structure 
(Fig.  1),  in  conflict  with  the  expectation  of  a  simple- 
rectangular  surface  crystal  structure  of  Sb  atoms.  In 
fact  the  atom  at  the  face-center  is  brighter  than  the 
other  four  atoms  of  the  unit  cell  at  some  bias  voltages, 
and  appears  to  be  an  above-surface  interstitial  atom  — 
with  there  being  one  “interstitial”  per  surface  unit  cell. 
Under  negative  bias,  only  Sb  atoms  are  expected  to  be 
imaged;  but  these  “interstitials”  lie  neither  at  the  posi¬ 
tions  of  surface  In  atoms  nor  of  surface  Sb  atoms,  based 
on  the  expected  “rigid  rotation”  geometry*  of  the  (110) 
surface  relaxation  of  InSb  [1-7]. 

By  increasing  the  magnitude  of  the  bias  voltage,  we 
are  able  to  reduce  the  strength  of  the  “interstitial”  im¬ 
age,  until  all  of  the  “interstitials”  disappear  simultane¬ 
ously  in  all  unit  cells.  Clearly  these  ethereal  images  of 
“interstitied”-like  features,  which  disappear  simultane¬ 
ously  in  every  unit  cell  as  the  magnitude  of  the  bias 
voltage  is  varied,  are  unlikely  to  be  images  of  actual  in¬ 
terstitial  atoms  on  the  surface. 

In  the  remainder  of  this  paper  we  describe  how  the 
data  are  tzdcen,  and  present  theoretical  arguments  to  the 
effect  that  the  “ethereal  interstitials”  are  in  fact  second- 
layer  Sb  atoms  that  become  rather  bright  in  comparison 
with  the  surface  Sb  atoms  (and  hence  appear  to  lie  above 
the  surface)  because  they  are  associated  with  states  near 
the  top  of  the  valence  band  of  InSb. 

2.  Experiment 

For  these  studies,  we  used  a  Pachyderm-4  cleav¬ 
age  scanning  tunneling  microscope,  which  has  been  de¬ 
scribed  in  detail  elsewhere  [Sj.  The  InSb  samples  were 
nominally  undoped  single-crystal  bars  manufactured  by 
Sumitomo.  Before  each  STM  experiment,  each  sample 


was  freshly  cleaved  in  UHV  (sslO~^®  torr),  producing  a 
mirror-like  flat  surface.  The  microscope  tips  were  made 
of  W  which  had  been  electrochemically  etched  in  an 
NaOH  solution. 

For  sample  biases  of  a  few  tenths  of  a  Volt,  the 
expected  rectangul2ir  surface  crystal  structure  was  ob¬ 
served  [8],  consisting  of  Sb  (In)  atoms  for  negative  (pos¬ 
itive)  sample  bias.  WTien  the  magnitude  of  the  negative 
bias  was  reduced  to  %0.05  V,  the  ethereal  interstitials 
appeared  at  the  center  of  each  stirface  unit  cell  (Fig. 
1).  The  intensity  of  this  feature  relative  to  the  inten¬ 
sity  of  the  Sb  corner-atom  intensity  v-aried  reproducibly 
from  zero  to  greater  than  unity  as  the  magnitude  of  the 
bias  voltage  was  reduced  from  a  few  tenths  of  a  Volt 
(or  more)  to  0.05  V.  The  smallest  and  best  bias  voltage 
magnitude  was  0.05  V';  below'  this,  no  features  could  be 
seen. 

3.  Theory 

The  ethereal  interstitial  features  lie  directly  above  the 
sites  of  second-layer  Sb  atoms  (Fig.  2),  and  so  we  inves¬ 
tigated  the  possibility  that  somehow  those  second-layer 
atoms  could  produce  an  image  that  appeared  to  have  an 
above- surface  interstitial  for  biases  of  order  -  0.05  V.  but 
would  have  that  feature  disappear  for  a  more-negative 
bias  of  -  0.3  V. 

To  do  this,  we  calculated  the  expected  images  for  bi¬ 
ases  of  -1.5  V  and  -  0.05  V,  producing  the  results  of  Fig. 
3.  Note  that  the  “ethereal  interstitial”  feature  does  ap¬ 
pear  prominently  in  Fig.  3b  (-  0.05  V),  but  not  in  Fig. 
3a  (-1.5  V),  and  therefore  is  explained  by  the  theory.  In 
fact.  Fig.  4  shows  the  intensity  ratio,  I2  of  the  “ethereal 
interstitial”  peak  to  Ij  for  the  1x1  peaks  (the  comers 
of  the  surface  unit  cell),  showing  that  the  theory  and 
the  data  exhibit  the  same  general  trends  as  the  bias  is 
varied. 

The  theoretical  calctilations  rely  on  B^deen’s  approx¬ 
imation  for  the  tunneling  current  [9]: 

I(V)  =  (2«/ft)  ^  f(E^)(l-f(E^-  cV)) 

(1) 

where  f(E)  is  the  Fermi  factor,  is  the  transfer  ma- 


trix  element  between  the  sample  state  v  and  the  tip 
state  and  V  is  the  bias  voltage.  Tersoff  and  Hamann 
(10]  have  obtained  an  approximate  expression  for  the 
tunneling  current,  \'alid  for  small  V: 

I(V)  «  V  2^  )|2<(E„-Ef), 

(2) 

where  tyi/(Rtip)  electronic  wavefunction  of  the 

sample  state  v  at  the  tip  position  R^jpi  and  Ep  is  the 
Fermi  energy.  The  summation  on  the  right  hand  side 
of  Eq.  (2)  is  basically  an  energj'-selected  charge  den¬ 
sity  (for  filled  states)  or  a  virtual  energj'-selected  charge 
density  (for  empty  states)  at  the  tip  position,  and  was 
obtained  by  Taylor  expansion  of  Bardeen’s  formula  to 
lowest  order  in  V.  As  such,  the  Tersoff-Hamann  expres¬ 
sion  allows  tunneling  for  states  at  Ep  only  and  does  not 
contain  the  Fermi  factors  of  Eq.  (1),  which  allow  tun¬ 
neling  for  states  with  energies  between  Ep  and  Ep  -  eV. 
Nevertheless  it  is  reasonable  to  extend  Eq.  (2)  to  higher 
bias  voltages  by  replacing  the  ^-function  with  a  window 
function  that  is  unity  within  the  interval  (Ep,Ep  -  eV) 
and  zero  elsewhere.  We  fxurther  simplify  our  calculation 
by  assuming  that  the  STM  image,  which  is  measured 
under  constant- current  conditions,  is  adequately  simu¬ 
lated  by  the  energj* -selected  charge  density  at  a  constant 
tip-height  of  5  A  above  the  surface. 

The  computations  were  executed  using  the  local- 
density  approximation  as  implemented  with  the  well- 
established  pseudofunction  method  (11)  —  a  method 
which  has  successfully  described  the  electronic  and  op¬ 
tica)  properties  of  many  solids  [12]  and  solid  stirfaces 
(13].  Details  of  the  calculational  approach  are  available 
in  Ref.  (11]. 

4.  Results 

Because  of  the  p-like  anion-like  character  of  the  elec¬ 
tronically  filled  states  of  zinc- blende  semiconductors, 
STM  images  of  (110)  surfaces  for  negative  sample  bias 
(namely  for  tunneling  of  electrons  from  the  valence  band 
to  the  STM  tip)  are  expected  to  show  the  anion  atoms 
at  the  surface,  and  these  exhibit  a  rectangular  unit  cell 
for  the  (110)  surface.  Under  such  bias  conditions,  the 
cation  atoms  are  not  normally  visible.  Thus  the  charac- 


terislic  signature  of  tunneling  from  a  zinc- blende  (110) 
surface  is  a  rectangular  unit  cell  of  anion  atoms.  This 
signature  has  been  observed  in  the  STM  data  for  the 
InSb  (110)  surface  [8]  and  the  (110)  surfaces  of  other 
III-V  zinc-blende  semiconductors  (14),  for  moderate  to 
large  bizis  magnitudes,  and  in  the  theory  of  Fig.  3. 

The  striking  fact  is  that  the  theorj'  also  reproduces 
the  unexpected  “ethereal  interstitial”  feature  for  small 
negative  biases,  and  that  this  feature  is  associated  with  a 
directed  p-orbital  centered  on  a  sub-surface  Sb  atom  [15- 
17j.  In  the  early  days  of  scanixing  tunneb'ng  microscopy 
it  was  generally  believed  that  the  STM  image  is  a  topo¬ 
graph,  giving  a  picture  of  the  surface  atoms,  with  image 
intensities  reflecting  heights  indicative  of  the  actual  sur¬ 
face  atoms’  sizes.  This  belief  seems  to  be  held  by  some 
workers  even  today.  The  ethereal  interstitials  necessar¬ 
ily  fall  outside  this  notion  of  topographic  images,  be¬ 
cause  they  correspond  to  second-layer  Sb  atoms  that  are 
expected  to  be  invisible  due  to  their  lack  of  height,  and 
yet  are  observed  to  be  brighter  than  first-layer  atoms 
at  small  bias  voltages.  By  dissecting  our  calculations, 
we  find  that  the  relative  prominence  of  the  “ethereal  in¬ 
terstitials”  is  due  to  the  relative  amplitudes  of  surface 
and  sub- surface  Sb  wavefunctions  at  near-zero  bias.  The 
tunneling  conditions  for  small  negative  bi«  are  such 
that  only  electrons  from  near  the  top  of  the  InSb  va¬ 
lence  band  (at  k=0)  participate  in  the  tunneling.  At 
k=0  the  wavefunction  is  quasi-uniform  throughout  the 
crystal  and  one  would  expect  the  first-layer  w'avefunc- 
tion  to  produce  a  brighter  image,  were  it  not  for  the 
following:  (i)  The  surface  sp^  bonds  re-hybridize  as  a 
result  of  surface  relaxation  [16)  and  their  contribution 
to  the  tuimeling  in  the  energ\*  range  between  Ep  and 
Ep  -  0.05  eV  is  reduced  because  they  dangle  into  the 
vacuum  less  in  this  spectral  region.  In  contrast,  the 
more  distant  sub- surface  bonds  still  point  toward  the 
surface  and  contribute  strongly  to  the  turmeling;  and 
(ii)  The  surface  atoms’  waveftmctions  have  eN'anescent 
tails  which  decay  exponentially  away  from  the  surface, 
whereas  the  sub-surface  atoms  have  wavefunctions  with 
Wigner-Seitz  boundary-condition  dangling  tails  that  are 
not  appreciably  attenuated  until  they  reach  the  surface. 
Primarily  the  first  of  these  two  effects  causes  the  wave- 
function  amplitudes  of  the  surface  atoms  to  be  smaller, 


and  is  responsible  for  the  darkening  surface  atom  im¬ 
ages  and  the  appearance  of  the  second-layer  atoms  as 
relatively  brighter.  But  as  the  wavefunctions  of  the 
tunneling  electrons  become  less  k=0-like  on  the  aver¬ 
age,  namely,  as  the  bias  voltage  becomes  more  negative, 
the  dangling-bond  states  at  the  surface,  which  are  spec- 
treilly  concentrated  at  energies  in  the  valence  band,  con¬ 
tribute  more  to  the  tunneling  and  eventuzdly  dominate 
it.  Then  the  “ethereal  interstitial”  becomes  (relatively) 
less  prominent,  with  its  relative  intensity  exhibiting  the 
same  trend  with  bias  voltage  as  the  theory  (Fig.  4).  Evi¬ 
dence  supporting  this  viewpoint,  that  the  surface  atoms 
are  darker  in  the  image  due  mostly  to  their  being  re¬ 
hybridized  at  a  relaxed  surface,  is  that  the  ethereeil  in¬ 
terstitials  are  not  visible  in  the  image  computed  for  the 
unrelaxed  surface  under  the  same  bias  conditions. 

Frequently  one  expects  a  topographic  image  from 
an  STM  experiment  [18],  but  non-topographic  images 
with  bright  second-layer  atoms,  similar  to  what  we  have 
found  here,  do  occur;  Hamers  et  al.  [19]  have  re¬ 
ported  images  with  sub-surface  features  for  the  clean 
Si  (111)  7x7  surface.  Tsai  ei  al.  have  reported  j5-SiC 
(111)  surfaces  in  which  the  carbonaceous  surface  forms 
a  graphitic  monolayer,  and  the  STM  images  are  domi¬ 
nated  by  tunneling  from  those  second-layer  Si  dangling 
bonds  that  are  not  masked  bv  the  graphitic  surface  layer 
120). 


5.  Summary 

We  have  observed  InSb  (110)  surface  structures  by 
STM  under  negative  bias  that  feature  what  appear  to 
be  interstitial  atoms  in  every  unit  cell,  and  the  inter¬ 
stitials  appear  to  lie  above  the  surface  defined  by  the 
four  Sb  atoms  in  the  surface  unit  cell.  These  apparent 
interstitials  su'e  in  fact  sub-surface  Sb  atoms  that  are 
relatively  prominent  in  the  STM  images  taken  at  small 
negative  bias  voltages,  but  disappear  as  the  magnitude 
of  the  bias  is  increased. 

Our  results  illustrate  the  fact  that  non-topographic 
images  can  appear  —  and  can  be  very  deceptive  unless 
one  carefully  compares  theory  with  data. 
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FIGURE  CAPTIONS 


Fig.  1.  STM  image  of  the  InSb  (110)  surface  taken  at 
a  sample  bias  of  -  0.05  V.  Note  the  prominent  “inter¬ 
stitial”  at  the  face-center  of  the  unit  cell,  labeled  by  -f 
and  surroimded  by  less  bright  Sb  atoms  on  the  corners 
of  a  rectangular  tmit  cell. 


Fig.  2.  Top  view,  looking  down,  of  the  (110)  InSb  sur¬ 
face,  with  the  Sb  atoms  (solid  circles)  and  the  In  atoms 
(open  circles).  The  surface  unit  cell  is  dashed.  The 
second-layer  Sb  is  denoted  by  a  small  solid  circle,  and 
is  directly  below  the  apparent  position  of  the  “ethereal 
interstitial.” 


Fig.  3.  Calculated  STM  image  of  the  InSb  (110)  surface,  as  discussed  in  the 
text,  (a)  Image  for  a  sample  bias  of -1.5  V,  showing  the  expected  1x1  rectangular 
unit  cell,  with  an  Sb  atom  slightly  to  the  left  of  the  “x”.  (b)  Image  for  a  bias 
of  -  0.05  V,  showing  the  “ethereal  interstitial”  at  the  face-centered  site,  with  the 
“interstitial”  just  to  the  left  of  the  “-I-”  and  a  comer  atom  of  the  1x1  cell  to  the 
left  of  the  “x”. 


Fig.  4.  Theoretical  (circles)  and  experimental 
(squares)  intensity  ratios  I2/I1  versus  bias  voltage  (in 
V)  for  the  InSb  (110)  surface.  I2  is  the  intensity  of  the 
“ethereal  interstitial”  peak,  and  Ij  is  the  intensity  of  the 
corner  Sb  atoms  of  the  1x1  surface  unit  cell. 
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ABSTRACT 

Three  types  of  steps  are  observed  on  the  cleaved  InP(llO)  surf^ace,  using  atomic- 
resolution  ultra-high  vacuum  (UHV)  scanning  tunneling  microscopy  (STM).  The  step 
edges  are  oriented  along  the  (110),  (111),  and  (112)  directions.  Atomic  models  of 
monatomic- height  (111)  and  (112)  steps  indicate  that  the  edges  of  each  of  these  unre¬ 
lated  steps  should  have  pairs  of  dangling  bonds.  We  propose  that  the  bonds  dimerize, 
causing  the  edges  to  relax  and  form  periodic  structures  along  the  edge. 

I.  INTRODUCTION 

In  this  paper  we  present  scanning  tunneling  naicroscopy  (STM)  images  of  the 
InP(llO)  siirface  cleaved  in  ultra- high  vacuum,  and  find  atomic-resolution  images  of 
steps  on  this  surface.  Ebert  et  al  [1]  have  examined  the  InP(llO)  surface  with  STM, 
but  the  present  paper  is,  we  believe,  one  of  the  first  STM  efforts  (2]  to  study  such 
steps,  which  are  important  as  nucleation  sites  for  layered  InP(llO)  growth,  very  likely 
determine  essential  physics  of  contact  metallization,  doubtless  play  a  role  in  Schottky 
barrier  formation,  and  getter  defects.  Here  we  present  images  of  the  flat  (110)  surface 
and  of  steps  on  the  surface.  We  also  propose  a  Dimerization  Model  of  atomic  relaxation 
at  (112)  step  edges  on  the  (110)  siirface. 

II.  EXPERIMENT 

The  images  were  obtained  with  a  Pachyderm-4  ultra-high  \'acuum  (UHV)  scanning 
tunneling  microscope  equipped  with  a  demountable  sample  holder  containing  six  sam¬ 
ples  at  one  time.  Each  sample  was  cleaved  in  UHV  (10"^^  torr  range)  before  being 
imaged  with  the  tungsten  tip  of  the  STM.  The  microscope,  cleavage  method,  and  tip 
preparation  are  described  in  detail  elsewhere  [3]. 

III.  RESULTS  AND  DISCUSSION 

We  imaged  the  (110)  surface  of  InP  with  atomic  resolution.  An  example  of  aflat 
(110)  surface  imaged  by  tunneling  into  the  empty  conduction  band  states  is  shown  in 
Fig.  1.  Similar  images  for  the  InSb(llO)  surface  may  be  found  in  Ref.  (3).  In  all 
cases  the  expected  1x1  rectangular  lattice  was  observed,  and  only  one  atomic  species 
was  imaged  at  a  time,  depending  on  the  polarity  of  the  sample  bias*,  indium  atoms  for 
positive  bias,  and  phosphorus  atoms  for  negative  bias. 

In  addition  to  flat  terraces,  steps  on  these  surfaces  were  also  observed.  For  InP(llO), 
steps  were  found  with  edges  running  in  the  (110),  (HI),  and  (112)  directions.  An  ex¬ 
ample  showing  sJl  three  step  types  on  InP(llO)  is  shown  in  Fig.  2,  where  a  single 
triangular-shaped  terrace  was  formed  by  intersecting  (110),  (HI),  and  (112)  steps. 
This  intersection  of  steps  on  the  (110)  surface  of  InP  was  unanticipated.  These  steps 
are  the  result  of  the  cleavage  process  and  do  not  represent  an  equilibrium  step  con¬ 
figuration,  as  for  some  heat-treated  samples.  On  Si  and  Ge  (001)-2xl  and  (lll)-7x7 
surfaces,  the  steps  are  normally  created  not  by  cleavage  but  by  high  temperature  pro¬ 
cessing.  namely  under  quasi-equilibrium  conditions.  Our  steps  on  the  InP(llO)  surface 


Fig.  1,  STM  image  showing  empty  states  on  the  flat,  cleaved  InP(llO)  surface:  A 
35  Ax 35  A  image  taken  under  +2.0  V  sample  bias  with  a  tunneling  current  of  100  pA. 


Fig.  2.  Triangular  terrace  on  the  InP(llO)  surface,  observed  by  STM  at  a  voltage 
of  +2.0  V  and  a  current  of  100  pA.  A  triangular  structure  formed  by  intersecting 
(110),  (111),  and  (112)  steps  on  the  InP(llO)  surface  is  pictured.  The  region  shown  is 
450  Ax 450  A. 
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Fig.  3.  Dimerization  Model  for  a  (112)-type  step  edge,  the  (lT2)  step  edge,  before 
any  relaxation,  (a)  Top  view  and  (b)  side  view  of  the  (lT2]  step  showing  dangling  bonds 
(“rabbit  ears”)  before  dimerization. 
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Fig.  4.  Dimerization  Model  for  a  (112)-type  step  edge,  the  (112)  step  edge,  after 
dimeriiation.  For  simplicity  of  presentation,  the  29®  rigid  rotation  (6)  of  the  P  anions 
up  out  of  the  surface  has  been  omitted,  (a)  Top  view  and  (b)  side  view  of  the  [lT2) 
step  showing  dangling  bonds  after  dimerization. 


are  not  created  in  quasi-equilibrium,  but  by  cleavage  (a  rather  non-equilibrium  event). 
If,  on  the  one  hand,  step-intersections  on  the  InP(llO)  surface  were  to  occur  at  equilib¬ 
rium,  then  those  steps  would  have  effective  attractive  interactions.  On  the  other  hand, 
if  the  step  intersections  are  a  result  of  cleaving,  then  the  physics  of  their  intersections 
should  produce  insights  into  the  processes  by  which  the  cleaving  energy  is  dispersed. 

Simple  geometric  models  of  one-atom-high  steps  can  be  constructed  from  the  zinc- 
blende  (110)  surface  structure.  Fig.  3  shows  a  (ll2)-type  step,  assumed  to  be  one-atom 
in  height,  both  top  and  side  views.  Note  that  there  are  two  dangling  bonds  periodically 
repeated  at  the  step-edge,  one  parallel  to  the  surface  and  another  directed  at  an  angle 
of  55®  with  respect  to  the  surface  plane. 

We  propose  that  these  bonds  will  dimerize  and  the  stepj-edge  will  relax,  forming 
the  surface  structure  of  Fig.  4.  One  implication  of  this  Dimer  Model  of  the  (ll2)-t>'pe 
monatonnic  step  is  that  the  relaxed  step-edge  should  be  scalloped,  showing  a  periodic 
structure  of  about  7.2  A.  We  also  propose  that  the  (111)  monatomic  step  exhibits  sim¬ 
ilar  dimerization,  which  should  lead  to  a  periodicity  of  10.2  A  on  the  InP(l  10)  surface. 
At  the  present  we  are  unable  to  determine  whether  the  dimerized  atoms  in  monatomic 
steps  are  In  pairs  or  P  pairs;  theoretical  work  on  this  question  is  in  progress. 

In  recent  studies  of  the  InAs(llO)  surface,  we  have  searched  for  and  found  such 
periodicities  on  monatomic  steps  —  which  will  be  reported  soon  [4).  Presumably  the 
seune  edge-reconstruction  effect  will  be  present  in  InP.  We  note  also  that  not  all  steps 
observed  are  monatomic,  and  that  models  of  these  larger  steps  need  to  be  developed 

(5). 

It  is  worth  remarking  that  the  InP(llO)  surface  is  considerably  more  difficult  to 
image  with  STM  than  either  InAs(llO)  or  InSb(llO).  This  was  unanticipated,  and, 
while  at  first  we  thought  it  might  have  been  related  to  the  quality  of  our  samples,  we 
now  suspect  that  it  may  be  related  to  the  ionicity  of  the  material,  since  InSb  was  easy 
to  image  and  InAs  intermediate.  We  currently  have  no  explanation  of  this. 

IV.  CONCLUSIONS 

For  cleaved  InP(llO)  surfaces  we  observed  the  expected  1x1  sxirface  reconstruction 
which  is  consistent  with  the  standard  model  of  this  surface;  the  anions  rotate  up  out 
of  the  surface,  almost  in  a  rigid  rotation,  by  about  29®  [6).  In  addition  we  observed 
steps  running  in  preferred  directions  on  this  surface:  Step  edges  were  found  to  nm  in 
the  (110),  (111),  and  (112)  directions.  Model  calculations  axe  underway  to  determine 
the  energetics  of  formation  of  the  various  steps. 

Periodicity  along  the  (111)  and  (112)  steps  is  a  natural  consequence  of  the  Dimer 
Models  we  have  proposed  for  the  step-edges’  relaxation.  Model  calculations  are  un¬ 
derway  to  determine  if  the  dimers  along  the  (111)  and  (112)  steps  are  III-III  or  V-V 
dimers.  The  study  of  these  steps  appears  to  be  a  very  interesting  scientific  problem 
with  technological  implications. 
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Abstract 

Clean  surfaces  of  Si(100)-(2xl)  were  bombarded  by  3keV  Ar'*'  ions  at  doses  10**- 
10*^  ions  cm'^.  Scanning  tunneling  microscopy  (STM)  images  of  the  bombarded  surfaces  i 

showed  random  defects  in  the  form  of  missing  dimers  or  'craters*.  These  defects  ordered 
into  line  defects  perpendicular  to  the  dimerrows  upon  annealing  the  surface  at  temperatures 
between  600  °C  and  850  for  2  minutes.  Molecular  Dynamics  simulations  were  per- 

1 

formed  to  explain  the  shapes  and  sizes  of  the  observed  craters  and  also  to  examine  the  sta¬ 
bility  of  ordered  defects.  Our  calculations  showed  a  small  energy  difference  between  line 
defects  along  and  perpendicular  to  the  dimer  rows,  with  a  lower  energy  favouring  the  latter, 
in  agreement  with  STM  observations. 

PACS  numbers:  61.80.Jh;  68.55.Ln;  79.20.Rf  L 

i: 


1 


In  a  recent  experiment  in  which  we  bombarded  Si(100)-(2xl)  surfaces  with  3keV  Ar"*" 
ions  at  low  doses,  i.e.  -10^^  ions  cm'^  or  -10*^  monolayer  (ML),  at  room  temperature,  we 
reponed  the  observation  of  random  defects  which  ordered  into  line  defects  upon  annealing 
the  bombarded  surfaces  to  7S0°C  for  2  minutest  In  another  experiment  by  Bedrossian  and 
Klitsner^,  Si(100)'(2xl)  surfaces  were  bombarded  by  225  eV  Xe”*"  ions  at  high  doses,  i.e. 
lO^^  ions  cm‘^,  at  450  °C.  They  observed  retraction  of  the  Sg  edges,  and  subsequent 
appearance  of  vacancy  islands  on  the  A  terraces  upon  further  sputtering.  Here  we  use 
Chadi’s  notation^  of  and  Sg  refenring  to  step  edges  parallel  and  perpendicular  to  the 
dimer  rows  of  the  upper  terrace  respectively,  and  A  and  B  terraces  referring  to  those 
domains  with  dimer  rows  parallel  and  perpendicular  to  the  step  edges  respectively.  Our  low 
dose  experiment  is  illustrated  in  Fig.l.  The  sianing  Si(100)-(2xl)  surface  was  remarkably 
free  of  any  metal  contamination  which  could  give  rise  to  line  defects  perpendicular  to  the 
dimer  rows^.  Fig.  1(a)  shows  an  image  of  such  a  surface  taken  by  scanning  nmneling  micro¬ 
scopy  (STM),  with  single  missing  dimers  being  the  predominant  defects,  fig.  1(b)  shows 
the  surface  after  bombardment  by  3keV  Ar'*'  ions  at  normal  incidence  with  a  dose  of  1.5  x 
10  ions  cm*'^.  A  higher  concentration  of  random  defects  was  observed,  some  of  which 
contain  multiple  missing  dimers.  When  a  bombarded  surface  was  annealed  at  600  ®C  for  2 
minutes,  the  random  defects  appear  to  order  into  line  defects  perpendicular  to  the  dimer 
rows,  as  shown  in  Fig.  1(c).  Annealing  at  higher  temperatures,  750 ®C  and  850  °C  for  2  min¬ 
utes,  has  the  same  effect,as  shown  in  Fig.  1(d).  When  the  bombarded  surface  was  annealed 
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to  1200  °C  for  several  seconds,  the  surface  was  completely  restored  to  that  shown  in 
Fig.  1(a),  without  any  trace  of  line  defects.  We  did  not  observe  ordering  of  the  random 
defects  when  we  annealed  the  surface  at  500  for  2  minutes.  This  implies  that  the  thresh* 

old  temperature  for  ordering  lies  between  500  °C  and  600  ®C,  which  is  higher  than  the  450 
observed  by  Bedrossian  and  Klitsnei^  for  vacancy  islands  formation  during  high  dose 
bombardment.  This  also  implies  that  we  have  studied  the  formation  and  ordering  of  random 
defects  separately.  Interestingly,  the  vacancy  islands,  which  have  much  larger  areas  than  the 
line  defects,  have  their  long  edges  parallel  to  the  dimer  rows,  in  contrast  to  our  line  defects 
which  are  perpendicular  to  the  dimer  rows. 

In  this  report,  we  compare  our  STM  images  with  the  results  of  molecular  dynamics 
(MD)  simulations  in  order  to  understand  the  mechanisms  which  lead  to  the  formation  of  the 
random  and  ordered  defects  on  the  surface  bombarded  by  low  dose  Ar'*’  ions. 

In  the  classical  MD  simulations  we  have  bombarded  a  dimerized  Si(lOO)  surface  with 
individual  3  keV  Ar'*’  ions.  The  microciystallite  consists  of  12  layers  of  32  atoms  each.  Peri* 
odic  boundary  conditions  are  applied  in  order  to  minimize  boundary  effects.  The  bottom 
layer  atoms  are  coupled  to  their  initial  positions  with  a  harmonic  force.  The  value  of  the 
force  constant  is  10'^  Nm'^  A  temperature  control  is  added  to  mimic  the  heat  flow  to  the 
bulk.  The  method  described  by  Berendsen  et  al^.  is  used  to  couple  all  the  atoms  of  the  crys¬ 
tallite  to  an  external  bath  at  constant  temperature.  In  the  simulation,  the  coupling  is  to  a  bath 
at  7  s  300^  with  a  time  constant  x  ss  400/s.  With  this  value  the  crystallite  is  cooled 
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down  to  300  K  within  a  few  ps  after  the  ion  impact.  In  order  to  coirect  for  the  dissipation 


of  energy  due  to  the  shock  wave  shortly  after  the  impact,  the  potential  energy  build  up  in 
the  bonom  atoms  is  removed  from  these  atoms  as  long  as  their  temperature  exceeds  the 
temperature  of  the  external  bath.  In  case  of  events  causing  large  surface  damage,  a  larger 
crystallite  is  used  containing  128  Si  atoms  in  each  layer.  This  size  is  large  enough  to  prevent 
interference  due  to  the  periodic  boundary  conditions  and  the  finite  size  of  the  crystallite. 
The  Si-Si  interactions  are  described  by  the  Stillinger  and  Weber  many-body  potential®.  To 
incorporate  a  more  repulsive  Si-Si  interaction  at  short  intemuclear  separations,  the  poten¬ 
tial  is  splined  between  r  *1.09  A  and  r  *  1.59  A  to  a  Moliere  potential  with  a  screening 
length  of  0.83  times  the  Firsov  value^.  The  attractive  Ai-Ar  and  Si- Ar  interactions  are  neg¬ 
ligible  under  the  experimental  conditions.  The  non-bonding  Ar-Si  and  Ar-Ar  interactions 
are  described  by  a  Moliere  potmdal,  also  with  a  screening  length  of  0.83  times  the  Firsov 
value. 

A  set  of  30  trajectories  is  developed  by  random  sampling  of  the  representative  area 
of  the  Si(lOO)  surface.  With  every  new  trajectory,  a  fresh  surface  is  bombarded  with  the  ion. 
This  approach  simulates  the  low  dose  experiment.  In  the  simulated  30  events  we  have 
observed  21  sputtered  atoms  and  57  Si  atoms  escaped  through  the  bottom  layer  of  the  crys¬ 
tallite.  The  sputter  yield  compares  well  with  the  MD  data  of  Smith  et  al.®  and  is  somewhat 
lower  than  the  experimental  high  dose  yield  of  ~1  atom/ion^.  This  implies  that  only  about 
one-third  of  the  surface  vacancies  are  created  by  sputter  removal.  The  remaining  vacancies 
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are  fonned  because  Si  atoms  in  the  top  layers  are  pushed  into  the  bulk. 

In  six  cases  the  simulation  produced  no  surface  damage.  In  the  remaining  simulations, 
we  observed  altogether  33  missing  dimers,  25  defect  dimers,  19  adatoms  and  1 1  adclusters. 
The  damage  is  in  the  form  of  17  small  crattrs,  8  medium  size  craters  and  1  big  crater.  In  a 
few  cases  the  ion  impact  created  more  than  one  crater  as  shown  in  Fig.2(a).  A  small  crater 
corresponds  to  the  area  of  abou  tone  dimer.  A  medium  crater  corresponds  to  a  damaged  area 
of  about  2  dimers  while  a  large  crater  has  an  area  of  about  8  dimers.  The  classification  is 
rather  arbitrary  and  sometimes  it  is  difficult  to  classify  a  crater  in  any  one  of  the  three  cat¬ 
egories.  Examples  of  small,  medium  and  large  craters  are  given  in  Fig.2(b.c,d). 

Our  simulations  produced  26  craters  from  30  ion  iizq)act5,  i.e.  ^1  crater  per  intact. 
From  our  STM  images  in  Fig.l,  the  clean  surface  in  Fig.l(a)  shows  -60  defects,  repre¬ 
sented  by  dark  spots  in  the  image;  while  the  bombarded  surface  in  Fig.l(b)  shows  -125 
defects  or  dark  spots  of  various  sizes.  There  is  an  increase  of  -65  defects  after  bombard- 
menu  The  ion  dose  in  Fig.l(b)  was  1.5  x  10^^  ions  cm'^,  which  translates  into  ~10  ion 
impacts  on  an  area  of  the  size  shown  in  Fig.l  (b).  So  the  number  of  additional  defects  is  six 
times  that  of  the  ion  impacts.  From  our  STM  images  taken  on  Si(l00)-(2xl)  surfaces  bom¬ 
barded  by  doses  varying  from  7xl0^^o6xl0*^  ions  cm’^,  we  observed  consistently  more 
defects  produced  by  the  bombardment  than  the  number  of  ion  impacts.  This  sharply  con¬ 
trasts  our  observations  on  the  Si(  1 1  l)-(7x7)  surfaced  on  which  we  found  the  ratio  of  defects 
(or  craters)  to  ion  impacts  was  -1.  The  origin  of  this  discrepancy  could  be  that  some  of  the 
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defects  produced  on  the  Si(100)*(2xl)  surface  may  be  electronic  in  origin  rather  than  top* 


ographic.  which  cannot  be  shown  in  our  classical  MD  simulations. 

It  is  interesting  to  compare  the  large  craters  observed  in  the  STM  images  with  the 
large  simulated  crater  shown  in  Fig.2(b).  Fig.3  shows  two  large  experimental  craters  on  the 
same  scale  as  that  in  Fig.  2(b).  There  is  a  remarkable  resemblance  between  experiment  and 
simulation.  In  the  simulated  crater,  we  observe  clusters  of  adatoms  in  the  vicinity  of  the  cra¬ 
ter.  We  also  observe  large  bright  spots,  which  may  be  due  to  adclusters,  near  the  large  cra¬ 
ters  in  the  STM  images. 

Annealing  at  temperatures  at  or  above  600  causes  the  random  defects  to  order  into 

narrow  line  defects  as  shown  in  Fig.  1(c)  and  (d).  The  origin  of  these  line  defects  may  be 
explained  if  we  assume  a  preferential  vacancy  mobility  along  the  dimer  rows^^,  similar  to 
the  mobility  of  Si  adatoms^  and  the  attractive  interaction  of  the  vacancies.  The  question 

of  whether  the  line  defects  should  be  perpendicular  or  parallel  to  the  dimer  tows  is  exam¬ 
ined  by  performing  a  microscopic  study  using  a  quantum-mechanical  MD  method.  The 
advantage  of  using  this  method  over  the  conventional  total  energy  calculadonal  method  is 
that  atomic  positions  are  automatically  optimized  according  to  Newton’s  laws  of  motion. 
In  contrast,  the  conventional  total  energy  methods  need  to  preconceive  some  prefen’ed 
geometries.  The  MD  method  employed  in  this  study  is  the  ab  iruiio  real-space  n>ethod  of 
Sankey  and  Niklewski^^,  which  has  been  shown  to  work  well  for  various  systems  of  carbon 
and  silicon*^'^^'*®. 
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In  this  study,  the  Si(lOO)  surface  is  modeled  by  a  five>Uyer>slab  repeated  slab  (super 
cell)  model;  the  atoms  in  the  bottom  layer  are  assigned  with  huge  masses  and  attached  with 
hydrogen-like  atoms  to  saturate  their  dangling  bonds  so  that  they  are  motionless  and  are 
regarded  as  bulk  atoms.  The  observed  bulk  latdce  constant  of  Si,  5.43072  A  is  used  to 
set  up  the  two-dimensional  lattice.  Chadi^^  has  shown  that  for  the  Si(lOO)  surface  atomic 
relaxations  extend  over  the  first  five  layers.  Since  the  atomic  displacement  for  the  fifth-layer 
atom  was  calculated  to  be  only  0.034  A  we  believe  the  use  of  a  fixed  fifth  layer  may  be 
adequate  for  studying  relative  preference  between  perpendicular  and  parallel  line  defects. 
(Here,  “perpendicular”  and  “parallel”  arc  defined  with  respect  to  the  dimer  rows).  The  line 
defect  is  modeled  by  a  periodic  array  of  infinitely  long  missing  atomic  lines.  Rg.  4(a)  and 
(b)  correspond  to  line  defects  which  are  perpendicular  and  parallel  to  dimo’  rows  respec¬ 
tively.  In  the  former  case,  denoted  as  case  (i),  the  missing  atomic  lines  are  seperated 
three  surface  atoms  (  or  three  pairs  of  dimers).  The  unit  cell  contains  thirty  eight  atoms,  i.e. 
six  in  the  surface  layer  and  eight  per  layer  in  the  second  to  fifth  layers.  As  shown  in  Fig.4(a), 
exposed  second  layer  atoms  have  freedom  to  form  dimers;  whether  they  form  dimers  is  to 
be  determined  by  the  MD  simulations,  which  let  all  atoms  except  those  in  the  bottom  layer 
settle  to  their  most  stable  positions.  As  for  the  case  shown  in  Fig.4(b),  denoted  as  case  (ii), 
the  twin  missing  atomic  lines  (or  missing  dimer  rows)  are  separated  by  three  dimer  rows. 
The  choice  of  twin  missing  atomic  lines  is  to  keep  the  ratio  between  the  number  of  missing 
atoms  and  surface  atoms  in  a  unit  cell  the  same  as  that  for  case  (i)  so  that  the  comparison 
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of  total  enei{ies  is  meaningful.  We  choose  a  unit  cell  twice  as  large  as  that  of  case  (i).  There 
axe  seventy-six  atoms  per  unit  cell;  twelve  in  the  surface  layer  and  sixteen  per  layer  in  the 
second  to  fifth  layers.  This  choice  of  a  doubled  unit  cell  is  to  allow  the  exposed  second- 
layer  atoms  freedom  to  form  dimers.  For  both  cases  (i)  and  (ii)  we  use  the  single  two- 
dimensional  special  point  of  Cunningham^^  for  a  rectangular  lattice  and  its  refiectionally 
symxoetric  points  with  respect  to  x  and  y  axes  and  the  inversion  symmetric  point  to  opfi- 
mize  atomic  positions.  To  compare  total  energies  (per  atom)  between  cases  (i)  and  (ii),  we 
find  that  the  use  of  sixteen  two-dimensional  special  points  of  Cunningham^^  is  quite  ade¬ 
quate.  We  have  also  tested  the  use  of  sixty-four  special  points  and  it  essentially  yields  the 
same  results. 

Fig.4(a)  and  (b)  show  the  lateral  displacements  of  surface-  and  second-layer  atoms 
from  their  ideal  unrelaxed  atomic  positions.  In  both  cases  (i)  and  (ii)  dimerization  directions 
of  surface  atoms  and  exposed  second-layer  atoms  can  be  seen  perpendicular  to  each  other 
as  stated  previously.  Surface  atoms  fonn  buckled  (not  shown  in  Fig.4)  asymmetric  dimers 
simular  to  those  discussed  by  Chadi^^.  In  case  (i).  the  second  layer  dimerization  is  not  very 
pronounced  and  the  atomic  pairs  have  bond  lenghts  (2.86A)  much  larger  than  that  of  sur¬ 
face  dimers  (2.36A).  This  is  because  each  exposed  second-layer  atom  is  still  bonded  to  and 
drawn  towards  a  surface  atom  at  either  edge  of  the  line  defect.  We  have  performed  charge 
density  calculations  and  find  that  the  charge  density  is  significantly  large  in  the  central 
region  of  each  atom-pair  in  the  exposed  layer  to  indicate  that  the  atom-pairs  indeed  form 
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dimers  through  their  single  dangling  bonds.It  is  interesting  to  note  that  in  case  (ii),  the  two 
second-layer  dimers  lines  in  the  defect  are  not  in  phase.  They  are  approximately  out  of 
phase  by  about  180^.  This  arrangement  exerts  less  strain  on  third-layer  atoms  bridged 
between  these  two  second-layer  dimer  lines  than  when  they  are  in  phase. 

We  find  that  the  total  energy  for  case  (i)  is  about  0.01 1  eV/atom  lower  than  that  for 
case  (ii)  favoring  line  defects  perpendicular  to  dimer  rows  as  shown  in  the  STM  images 
Fig.  1(c)  and  (d).  The  lower  total  energy  for  case  (i)  may  be  qualitatively  understood  from 
the  shorter  bonding  distance  ( and  hence  stronger  bonding)  between  an  exposed  second- 
layer  atom  and  a  surface  atom  at  either  edge  of  the  line  defect  as  shown  in  Fig.4(a).  How¬ 
ever,  the  preference  for  the  case  (i)  is  partly  compensated  by  the  restriction  of  exposed-sec- 
ond-layer  dimerization,  as  stated  earlier,  and  by  the  interruption  of  the  (2x1)  dimer  rows.  In 
contrast,  the  exposed  second-layer  atoms  in  case  (ii)  have  less  constraint  to  form  dimers  and 
have  no  interruption  of  dimer  rows  as  shown  in  Fig.4(b).  This  may  be  the  reason  why  the 
calculated  total  energy  difference  is  small. 

The  observations  by  Bedrossian  and  Klitsner^  appear  to  contradict  the  above  results 
of  our  calculations,  since  they  found  monolayer  deep  vacancy  islands  elongated  parallel  to 
the  dimer  rows.  However,  it  must  be  remembered  that  their  experiment  was  conducted  at 
450°C  with  Xe'*'  ions  bombarding  the  surface  at  a  dose  several  orders  of  magnitude  higher 
than  in  our  experiment.  In  their  case,  due  to  the  high  dose  bombardment,  the  mobility  of  the 
vacancies  is  enhanced  similar  to  the  enhanced  mobility  of  surface  atoms  in  ion-beam 
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assisted  deposition^.  Moreover,  the  mean  vacancy-migration  barrier  is  lower  due  to  the 
high  density  of  vacancies  with  fewer  neighbouring  dimers.  Both  effects  lead  to  the  estab- 
lisment  of  an  equilibrium  situation  without  invoking  any  presumpdon  of  anisotropic 
vacancy  diffusion .  The  elongation  of  their  vacancy  islands  along  the  dimer  rows  is  expected 
since  the  energy  of  the  step  is  lower  than  that  of  the  Sg  step^'*^. 

In  summary,  we  have  obtained  good  agreement  between  molecular-dynamics  simu- 
ladons  and  STM  images  of  random  and  ordered  defects  produced  by  eneigedc  ion  impacts 
on  the  Si(100)-(2xl)  surface. 

LS.T.  Tsong  thanks  Philips  Research  Laboratories  for  the  award  of  a  van  Houten  Fel¬ 
lowship.  M.-H  Tsai  and  J.D.  Dow  acknowledge  suppon  from  U.S.  Army  Research  Office 
( DAAL03-91-G-0054). 
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Figure  Captions 


Fig.  I  Filled-state  STM  images  of  Si(100)'(2xl)  surfaces  taken  at  sample  bias  •2V 
and  tunnel  cuirent  0.5  nA.  (a)  Oean  surface  piior  to  bonbaidment.  Scan  area 
s  28nm  x  28nm.  (b)  After  bombardment  by  3keV  Ar'*'  ions  at  a  dose  of  1.5  x 
10^^  ions  cm*^.  Scan  area  *  26nm  x  26nm.  (c)  Bombarded  surface  after 
annealling  for  2  minutes.  Scan  area  «  22nm  x  22nm.  (d)  Bombarded  surface 
after  annealing  at  750°C  for  2  minutes.  Scan  area  «  40nm  x  40nm. 

Fig.  2  Ion  impact  craters  produced  by  molecular  dynamics  simulations,  (a)  A  double 
crater,  (b)  A  large  crater,  (c)  A  medium  crater,  (d)  A  small  crater.  The  waters 
in  (b)  and  (c)  appear  double  because  of  periodic  boundary  conditions. 

Fig.  3  STM  images  of  two  large  craters  caused  by  a  3keV  Ar'*’  ion  bombardment. 
Area  »  6nm  x  6nm. 

Fig.  4  Molecular*dynamics  simulations  of  a  periodic  array  of  line  defects  on  the 
Si(l00H2xl)  surface,  (a)  Line  defect  perpendicular  to  dimer  rows,  and  (b) 
line  defect  parallel  to  dimer  rows.  The  small  crosses  and  the  solid  dots  are 
ideal  unrelaxed  positions  for  surface  and  second-layer  atoms  respectively,  the 
rectangle  enclosed  by  dashed  lines  is  the  unit  cell. 
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(Received  29  April  1992) 

Scanning  tunneling  microscopy  (STM)  images  of  Si(100)-(2x I)  surfaces  bombarded  by  low-dose  3- 
keV  Ar'*’  ions  showed  random  defects  which  ordered  into  line  defects  perpendicular  to  the  dimer  rows 
upon  annealing  at  elevated  temperatures.  Molecular-dynamics  simulations  were  performed  to  explain 
the  shapes  and  sizes  of  the  observed  random  defects  and  also  to  examine  the  stability  of  ordered  defects. 
Our  simulations  showed  good  agreement  with  STM  observations. 


PACS  numbers:  6l.80.Jh,  68  55.Ln.  79.20.Rf 

The  process  of  sputtering,  which  involves  the  removal 
of  individual  atoms  and  molecules  from  a  surface  by  en¬ 
ergetic  ions,  has  been  a  subject  of  interest  for  many  de¬ 
cades.  With  the  advent  of  scanning  tunneling  microscopy 
(STM),  one  has  the  means  of  viewing  the  surface  damage 
due  to  sputtering  with  atomic  resolution.  In  a  recent  ex¬ 
periment  in  which  we  bombarded  Si(l00)-(2x  i)  surfaces 
with  3-keV  Ar*^  ions  at  low  doses,  i.e.,  —10^^  ionsem 
or  — 10  monolayer  (ML),  at  room  temperature,  we  re¬ 
ported  the  observation  of  random  defects  which  ordered 
into  line  defects  upon  annealing  the  bombarded  surfaces 
to  750 for  2  min  II).  The  low-dose  experiment  is  illus¬ 
trated  in  Fig.  1.  The  starting  Si(l00)-(2x  |)  surface  was 
remarkably  free  of  any  metal  contamination  which  could 


FIG.  1.  Filled-state  STM  images  of  Si(100)-(2x  1)  surfaces 
taken  at  sample  bias  of  —2  V  and  tunnel  current  0.5  nA.  (a) 
Clean  surface  prior  to  bombardment.  Scan  area  is  28  nmX28 
nm.  (b)  After  bombardment  by  3-keV  Ar'*'  ions  at  a  dose  of 
1.5x10^^  ionsem'^  Scan  area  is  26  nmx26  nm.  (c)  Bom¬ 
barded  surface  after  annealing  at  600  for  2  min.  Scan  area 
is  22  nmx22  nm.  (d)  Bombarded  surface  after  annealing  at 
750 for  2  min.  Scan  area  is  40  nm  X40  nm. 


give  rise  to  line  defects  perpendicular  to  the  dimer  rows 
l2].  Figure  1  (a)  shows  an  STM  image  of  such  a  surface, 
with  single  missing  dimers  being  the  predominant  defects. 
Figure  Kb)  shows  the  surface  after  bombardment  by  3- 
keV  Ar*^  ions  at  normal  incidence  with  a  dose  of  1.5 
X 10*^  ionsem  equivalent  to  '—0.002  ML  coverage.  A 
higher  concentration  of  random  defects  was  observed, 
some  of  which  contain  multiple  missing  dimers.  When  a 
bombarded  surface  was  annealed  at  600  for  2  min,  the 
random  defects  appeared  to  order  into  line  defects  per- 
pcndicular  to  the  dimer  rows,  as  shown  in  Fig.  l(c).  An¬ 
nealing  at  higher  temperatures,  750  and  850 for  2 
min,  produced  the  same  effect  as  shown  in  Fig.  1(d). 
When  the  bombarded  surface  was  annealed  at  950 for 
2  min,  or  flashed  to  1200^C  for  a  few  seconds,  the  sur¬ 
face  was  completely  restored  to  that  shown  in  Fig.  I  (a), 
without  any  trace  of  line  defects,  thus  ruling  out  the  pos¬ 
sibility  of  metal  contamination  which  invariably  worsens 
with  each  annealing  cycle.  We  did  not  observe  any  order¬ 
ing  of  the  random  defects  when  we  annealed  the  surface 
at  500 for  2  min.  This  implies  that  the  threshold  tem¬ 
perature  for  ordering  lies  between  500  and  600 ®C.  Over 
100  STM  images  were  taken  in  our  experiment  and  con¬ 
sistent  results  were  obtained. 

The  possibility  exists  that  the  line  defects  are  not  due 
to  ordering  of  surface  vacancies,  but  arise  from  the  order¬ 
ing  of  implanted  Ar  atoms  diffusing  to  the  surface  during 
annealing.  However,  previous  ellipsometry  studies  of 
desorption  of  implanted  Ar  from  Si  surfaces  [3]  show 
that  desorption  takes  place  at  temperatures  as  low  as 
400 ^C,  reaching  a  maximum  at  600 ®C.  Complete  de¬ 
sorption  is  achieved  at  800 ‘^C  and  the  Si  surface  returns 
to  its  original  state  [4].  These  results  do  not  agree  with 
our  observation  that  ordering  occurs  between  600  and 
850 ®C,  and  the  line  defects  only  disappear  at  950 °C  or 
upon  flashing  the  surface  to  1200*^C.  Analysis  of  four¬ 
teen  STM  images  taken  before  and  after  annealing  shows 
that  the  number  of  vacancies  increases  on  average  by 
—  15%  after  annealing.'  A  plausible  cause  for  the  in¬ 
crease  is  the  migration  of  buried  vacancies  to  the  surface 
during  the  annealing  process,  or  it  could  be  the  result  of 
random  sampling  of  the  surface  because  the  STM  did  not 
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image  the  same  localized  area  after  annealing.  The 
Gaussian  distribution  of  the  Ar"^  ions  in  the  beam  spot 
impacting  the  surface  could  also  give  rise  to  the  counting 
discrepancy. 

In  this  Letter,  we  compare  our  STM  images  with  the 
results  of  molecular-dynamics  (MD)  simulations  in  order 
to  understand  the  mechanisms  which  lead  to  the  forma¬ 
tion  of  the  random  and  ordered  defects  on  the  surface 
bombarded  by  a  low  dose  of  Ar'*'  ions. 

In  the  classical  MD  simulations  we  have  bombarded  a 
dimerized  Si(lOO)  surface  with  individual  3-keV  Ar"^ 
ions.  The  microcrystallite  consists  of  12  layers  of  32 
atoms  each.  Periodic  boundary  conditions  are  applied  in 
order  to  minimize  boundary  effects.  The  bottom-layer 
atoms  are  coupled  to  their  initial  positions  with  a  har¬ 
monic  force.  The  value  of  the  force  constant  is  10”^ 
Nm”^  A  temperature  control  is  added  to  mimic  the 
heat  flow  to  the  bulk.  The  method  described  by  Ber- 
endsen  et  ai  [5]  is  used  to  couple  all  the  atoms  of  the 
crystallite  to  an  external  bath  at  constant  temperature. 
In  the  simulation,  the  coupling  is  to  a  bath  at  r**300  K 
with  a  time  constant  r  —400  fs.  With  this  value  the  crys¬ 
tallite  is  cooled  down  to  30  K  within  a  few  ps  after  the  ion 
impact.  In  order  to  correct  for  the  dissipation  of  energy 
due  to  the  shock  wave  shortly  after  the  impact,  the  poten¬ 
tial  energy  buildup  in  the  bottom  atoms  is  removed  from 
these  atoms  as  long  as  their  temperature  exceeds  the  tem¬ 
perature  of  the  external  bath.  In  the  case  of  events  caus¬ 
ing  large  surface  damage,  a  larger  crystallite  is  used  con¬ 
taining  128  Si  atoms  in  each  layer.  This  size  is  large 
enough  to  prevent  interference  due  to  the  periodic  bound¬ 
ary  conditions  and  the  finite  size  of  the  crystallite.  The 
Si-Si  interactions  are  described  by  the  Stillinger  and 
Weber  many-body  potential  [6l.  To  incorporate  a  more 
repulsive  Si-Si  interaction  at  short  internuclear  separa¬ 
tions,  the  potential  issplined  between  r“1.09  and  L59  A 
with  a  Moliere  potential  with  a  screening  length  of  0.83 
times  the  Firsov  value  [7l.  The  attractive  Ar-Ar  and  Si- 
Ar  interactions  are  negligible  under  the  experimental 
conditions.  The  nonbonding  Ar-Si  and  Ar-Ar  interac¬ 
tions  are  described  by  a  Moliere  potential,  also  with  a 
screening  length  of  0.83  times  the  Firsov  value.  The  indi¬ 
vidual  trajectories  are  sensitive  to  the  details  of  the 
short-range  interactions  such  as  small  variations  of  the 
impact  points  of  the  incoming  ion.  Consequently,  we  per¬ 
formed  statistical  analysis  on  the  simulations  in  much  the 
same  manner  as  we  did  on  the  experimental  results. 

A  set  of  30  trajectories  is  developed  by  random  sam¬ 
pling  of  the  representative  area  of  the  Si(lOO)  surface. 
With  every  new  trajectory,  a  fresh  surface  is  bombarded 
with  the  ion.  This  approach  simulates  the  low-dose  ex¬ 
periment.  In  the  simulated  30  events  we  have  observed 
21  sputtered  atoms.  The  sputter  yield  compares  well  with 
the  MD  data  of  Smith,  Harrison,  and  Garrison  [8]  and  is 
somewhat  lower  than  the  experimental  high-dose  yield  of 
—  1  atom/ion  [91.  Only  about  one-third  of  the  surface 
vacancies  are  created  by  sputter  removal.  The  remaining 


vacancies  are  formed  because  Si  atoms  in  the  top  layers 
are  recoiled  into  the  bulk.  These  recoil  Si  atoms  are  the 
beginning  of  collision  cascades  and  in  the  30  events,  57  Si 
atoms  escaped  through  the  bottom  layer  of  the  crystallite. 

In  six  cases  the  simulation  produced  no  surface  dam¬ 
age.  In  the  remaining  simulations,  we  observed  altogeth¬ 
er  33  missing  dimers,  25  defect  dimers,  19  adatoms,  and 
11  adclustcrs.  The  damage  is  in  the  form  of  17  small 
craters,  8  medium  size  craters,  and  1  large  crater.  In  a 
few  cases,  the  ion  impact  created  more  than  1  crater,  as 
shown  in  Fig.  2(a).  A  small  crater  corresponds  to  the 
area  of  about  1  dimer.  A  medium  crater  corresponds  to  a 
damaged  area  of  about  2  dimers,  while  a  large  crater  has 
an  area  of  about  8  dimers.  The  classification  is  rather  ar¬ 
bitrary  and  sometimes  it  is  difficult  to  classify  a  crater  in 
any  one  of  the  three  categories.  Examples  of  small,  medi¬ 
um,  and  large  craters  are  given  in  Figs.  2(b)- 2(d). 

Our  simulations  produced  26  craters  from  30  ion  im¬ 
pacts,  i.e.,  —  I  crater  per  impact.  From  our  STM  images 
in  Fig.  1,  the  clean  surface  in  Fig.  1(a)  shows  ^60  de¬ 
fects,  represented  by  dark  spots  in  the  image,  while  the 
bombarded  surface  in  Fig.  Kb)  shows  —125  defects  or 
dark  spots  of  various  sizes.  There  is  an  increase  of  —65 
defects  after  bombardment.  The  ion  dose  in  Fig.  Kb) 
was  1.5 X 10^^  ionsem"^  which  translates  into  —10  ion 
impacts  on  an  area  of  the  size  shown  in  Fig.  1(b).  So  the 
number  of  additional  defects  is  6  times  that  of  the  ion  im¬ 
pacts.  From  our  STM  images  taken  on  Si(l00)-(2xl) 
surfaces  bombarded  by  doses  varying  from  7xio“  to 
6x  10*^  ionscm“^  we  observed  consistently  more  defects 
produced  by  the  bombardment  than  the  number  of  ion 
impacts.  This  sharply  contrasts  our  observations  on  the 
Sid  1  l)-(7x7)  surface  [ll,  on  which  we  found  the  ratio 


FIG.  2.  Ion-impact  craters  produced  by  molecular-dynamics 
simulations,  (a)  A  double  crater,  (b)  A  large  crater,  (c)  A 
medium  crater,  (d)  A  small  crater.  The  craters  in  (b)  and  (c) 
appear  double  because  of  periodic  boundary  conditions. 
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of  defects  (or  craters)  to  ion  impacts  was  ^1.  The  origin 
of  this  discrepancy  could  be  that  on  Si(lOO),  each  ion  im¬ 
pact  creates  more  than  1  site  from  which  sputtering 
occurs,  as  MD  simulations  show  in  Fig.  2(a).  Some  of 
the  defects  generated  may  also  be  electronic  in  origin 
rather  than  topographic,  and  they  are  unrelated  to 
sputtering* 

It  is  interesting  to  compare  the  large  craters  observed 
in  the  STM  images  with  the  large  simulated  crater  shown 
in  Fig.  2(h).  Figure  3  shows  two  large  experimental 
craters  on  the  same  scale  as  that  in  Fig.  2(b).  There  is  a 
remarkable  resemblance  between  experiment  and  simula¬ 
tion.  In  the  simulated  crater,  we  observe  clusters  of  ada¬ 
toms  in  the  vicinity  of  the  crater.  We  also  observe  large 
bright  spots,  which  may  be  due  to  adclusters,  near  the 
large  craters  in  the  STM  images. 

Annealing  at  temperatures  at  or  above  600° C  causes 
the  random  defects  to  order  into  narrow  line  defects  as 
shown  in  Figs*  1  (c)  and  1  (d)*  A  possible  origin  of  these 
line  defects  could  be  a  preferential  mobility  of  vacancies 
on  the  surface  along  the  dimer  rows  [lOl  as  a  conse¬ 
quence  of  the  annealing  temperature,  and  the  attractive 
interaction  of  the  vacancies*  The  question  of  whether  the 
line  defects  should  be  perpendicular  or  parallel  to  the  di¬ 
mer  rows  is  examined  by  performing  a  microscopic  study 
using  a  quantum-mechanical  MD  method*  The  advan¬ 
tage  of  using  this  method  over  the  conventional  total- 
energy  calculational  method  is  that  atomic  positions  are 
automatically  optimized  according  to  Newton’s  laws  of 
motion*  In  contrast,  the  conventional  total-energy  meth¬ 
ods  need  to  preconceive  some  preferred  geometries*  The 
MD  method  employed  in  this  study  is  the  ab  initio  real- 
space  method  of  Sankey  and  Niklewski  [11],  which  has 
been  shown  to  work  well  for  various  systems  of  carbon 
and  silicon  [12-14]* 

In  this  study,  the  Si(lOO)  surface  is  modeled  by  a  five- 
layer-slab  repeated  slab  (supercell)  model;  the  atoms  in 
the  bottom  layer  are  assigned  with  huge  masses  and  at¬ 
tached  with  hydrogenlike  atoms  to  saturate  their  dangling 
bonds. so  that  they  are  motionless  and  are  regarded  as 
bulk  atoms*  The  observed  bulk  lattice  constant  of  Si, 
5*43072  A  [15],  is  used  to  set  up  the  two-dimensional  lat¬ 
tice*  Chadi  [16]  has  shown  that  for  the  Si(lOO)  surface 
atomic  relaxations  extend  over  the  first  five  layers.  Since 
the  atomic  displacement  for  the  fifth-layer  atom  was  cal¬ 


culated  to  be  only  0.034  A  [16],  we  believe  the  use  of  a 
fixed  fifth  layer  may  be  adequate  for  studying  relative 
preference  between  perpendicular  and  parallel  line  de¬ 
fects.  (Here,  ‘‘perpendicular”  and  “parallel”  are  defined 
with  respect  to  the  dimer  rows.)  The  line  defect  is 
modeled  by  a  periodic  array  of  infinitely  long  missing 
atomic  lines.  Figures  4(a)  and  4(b)  correspond  to  line 
defects  which  are  perpendicular  and  parallel  to  the  dimer 
rows,  respectively.  Because  of  the  usual  computer  con¬ 
straints,  we  could  realistically  only  consider  the  lower 
limit  of  the  line-defect  width.  In  the  perpendicular  case 
(i),  the  line  defect  is  a  single  missing  atomic  line  repeat¬ 
ing  itself  periodically  every  fourth  line*  The  unit  cell  con¬ 
tains  38  atoms,  i.e.,  6  in  the  surface  layer  and  8  per  layer 
in  the  second  to  fifth  layers.  As  shown  in  Fig.  4(a),  ex¬ 
posed  second-layer  atoms  have  the  freedom  to  form  di¬ 
mers;  whether  they  form  dimers  is  to  be  determined  by 
the  MD  simulations,  which  allow  all  atoms  except  those 
in  the  bottom  layer  to  settle  to  their  most  stable  positions* 
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FIG.  4.  Molecular-dynamics  simulations  of  a  periodic  array 
of  line  defects  on  the  Si(100)-(2x  J)  surface*  (a)  Line  defects 
perpendicular  to  dimer  rows,  case  (i);  and  (b)  line  defect  paral¬ 
lel  to  dimer  rows,  case  (ii).  The  small  crosses  and  the  solid  dots 
are  ideal  unrclaxed  positions  for  surface-  and  second-layer 
atoms,  respectively*  The  rectangle  enclosed  by  dashed  lines  is 
the  unit  cell. 
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For  the  parallel  case  shown  in  Fig.  4(b),  denoted  as  case 
(ii),  the  periodic  line  defects  are  in  the  form  of  a  single 
missing  dimer  row  (i.e.,  twin  missing  atomic  lines)  re¬ 
peating  every  fourth  dimer  row.  The  choice  of  twin  miss¬ 
ing  atomic  lines  is  to  keep  the  ratio  between  the  number 
of  missing  atoms  and  surface  atoms  in  a  unit  cell  the 
same  as  that  for  case  (i)  so  that  the  comparison  of  total 
energies  is  meaningful.  We  choose  a  unit  cell  twice  as 
large  as  that  of  case  (i).  There  are  76  atoms  per  unit  cell, 
12  in  the  surface  layer  and  16  per  layer  in  the  second  to 
fifth  layers.  This  choice  of  a  doubled  unit  cell  is  to  allow 
the  exposed  second-layer  atoms  freedom  to  form  dimers. 
For  both  cases  (i)  and  (ii)  we  use  the  single  two- 
dimensional  special  point  of  Cunningham  [l7]  for  a  rec¬ 
tangular  lattice  and  its  refiectionally  symmetric  points 
with  respect  to  x  and  y  axes  and  the  inversion  symmetric 
point  to  optimize  atomic  positions.  To  compare  total  en¬ 
ergies  (per  atom)  between  cases  (i)  and  (ii),  we  find  that 
the  use  of  16  two-dimensional  special  points  of  Cunning¬ 
ham  1 17]  is  quite  adequate.  We  have  also  tested  the  use 
of  64  special  points  and  it  essentially  yields  the  same  re¬ 
sults. 

Figures  4(a)  and  4(b)  show  the  lateral  displacements 
of  surface-  and  second-layer  atoms  from  their  ideal  unre¬ 
laxed  atomic  positions.  In  both  cases  (i)  and  (ii)  dimeri¬ 
zation  directions  of  surface  atoms  and  exposed  second- 
layer  atoms  can  be  seen  perpendicular  to  each  other  as 
stated  previously.  Surface  atoms  form  buckled  (not 
shown  in  Fig.  4)  asymmetric  dimers  similar  to  those  dis¬ 
cussed  by  Chadi  [I6l.  In  case  (i),  the  second-layer  di¬ 
merization  is  not  very  pronounced  and  the  atomic  pairs 
have  bond  lengths  (2.86  A)  much  larger  than  that  of  sur¬ 
face  dimers  (2.36  A).  This  is  because  each  exposed 
second-layer  atom  is  still  bonded  to  and  drawn  towards  a 
surface  atom  at  either  edge  of  the  line  defect.  We  have 
performed  charge-density  calculations  and  find  that  the 
charge  density  is  significantly  large  in  the  central  region 
of  each  atom  pair  in  the  exposed  layer  to  indicate  that  the 
atom  pairs  indeed  form  dimers  through  their  single  dan¬ 
gling  bonds.  It  is  interesting  to  note  that  in  case  (ii),  the 
two  second-layer  dimer  lines  in  the  defect  are  not  in 
ph  ase.  They  are  approximately  out  of  phase  by  about 
180°.  This  arrangement  exerts  less  strain  on  third-layer 
atoms  bridged  between  these  two  second-layer  dimer  lines 
than  when  they  are  in  phase. 

We  find  that  the  total  energy  for  case  (i)  is  about  0.01 1 
eV/atom  lower  than  that  for  case  (ii),  favoring  line  de¬ 
fects  perpendicular  to  dimer  rows  as  shown  in  the  STM 
images  Figs.  1(c)  and  1(d).  The  lower  total  energy  for 
case  (i)  may  be  qualitatively  understood  from  the  shorter 
bonding  distance  (and  hence  stronger  bonding)  between 
an  exposed  second-layer  atom  and  a  surface  atom  at  ei¬ 
ther  edge  of  the  line  defect  as  shown  in  Fig.  4(a).  How¬ 
ever,  the  preference  for  case  (i)  is  partly  compensated  by 
the  restriction  of  exposed-second-layer  dimerization,  as 
stated  earlier,  and  by  the  interruption  of  the  (2x  l)  dimer 
rows.  In  contrast,  the  exposed  second-layer  atoms  in  case 


(ii)  have  less  constraint  to  form  dimers  and  have  no  inter¬ 
ruption  of  dimer  rows  as  shown  in  Fig.  4(b).  This  may  be 
the  reason  that  the  calculated  total  energy  difference  is 
small. 

In  a  recent  report  by  Bedrossian  and  Klitsner  [18], 
Si(l00)-(2x  i)  surfaces  were  bombarded  by  225-eY  Xe"*^ 
ions  at  grazing  incidence  along  the  [llO]  direction  at 
temperatures  400-450  °C.  They  observed  retraction  of 
Sb  step  edges  as  well  as  monolayer-deep  depressions  or 
‘'vacancy  islands”  elongated  parallel  to  the  dimer  rows. 
However,  in  contrast  to  the  present  report,  the  structures 
they  observed  are  not  equilibrium  structures,  but  are 
structures  stabilized  at  moderately  high  temperatures, 
i.e.,  — 450°C,  by  kinetic  considerations,  analogous  with 
growth. 

In  summary,  we  have  obtained  good  agreement  be¬ 
tween  molecular-dynamics  simulations  and  STM  images 
of  random  and  ordered  defects  produced  by  energetic  ion 
impacts  on  the  Si(100)-(2x  l)  surface. 
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